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ABSTRACT 


The purpose of the Permanently Manned Autonomous Space Oasis, 
designated Project WISH: The Emerald City, is to serve as 
permanent living quarters for space colonists. In addition, it 
will serve as a stopover for space missions and will be capable of 
restationing itself practically anywhere within the solar system to 
provide support for these missions. The station should be self- 
sufficient, with no specific dependence on any resources from 
Earth. 

The 1990-1991 design team continued work started by last 
year’s class. This year, further studies were conducted in the 
areas of orbital mechanics, propulsion, attitude control, and human 
factors. Critical elements were identified in each of these areas, 
and guidelines were established for the design of the Emerald City. 
Using the knowledge gained from these studies, two particular 
missions of interest, a Saturn Envelope mission and an Earth-to- 
Mars mission, were examined. The size and mass estimates, along 
with the methodologies used in their determination, are considered 
to be the main accomplishments of Phase II. 
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FOREWORD 


This report reflects the result of our design classes through 
the second year of Project WISH. The design activity was carried 
out through a seguence of three courses offered during the autumn, 
winter, and spring quarters. The autumn quarter course was a 
technical elective, two credit-hour, introductory course, mainly 
dedicated to the understanding of the accomplishments of the first 
phase ( 1989 - 1990 ) of the project and building the background to be 
able to further the design. The winter quarter course fulfilled a 
four credit-hour design course requirement for the seniors. During 
the autumn and winter quarters, seminars by NASA Lewis personnel 
were included. The spring quarter class was a three credit-hour 
technical elective course which culminated in this report. Project 
WISH is measuring up to our expectations in every regard both in 
technical and educational aspects. 


Hayrani Oz 
Associate Professor 


iv 


TABLE OF CONTENTS 


INTRODUCTION x * * 

1.0 PROJECT OVERVIEW 1> 2 

1.1 PROJECT DESCRIPTION 1-3 

ORBITAL MECHANICS 2- 1 

2.0 INTRODUCTION 2.1 

2.1 NOMINAL ORBIT 2.2 

2.1.1 Background 2.2 

2.1.2 Mission Analysis 2.4 

2.1.3 Analysis Results 2.12 

Inner Planets 2.20 

Outer Planets 2.23 

2.1.4 Nominal Orbit Selection 2.24 

2.2 AV MINIMIZATION 2.27 

2.2.1 Introduction . . . • 2.27 

2.2.2 Background 2.28 

2.2.3 Assumptions 2.29 

2.2.4 Geometry 2.31 

2.2.5 Methodology 2.33 

2.3 CONCLUSIONS 2.36 

PROPULSION 3* 1 

3.0 INTRODUCTION 3.1 

3.1 GENERAL PROPULSION STUDY 3.2 

3.1.1 Theoretical Background 3.2 

3.1.2 General Propulsion Systems 3.5 

3.2 SPACE RADIATED GAS-CORE NUCLEAR ROCKET STUDY . . 3.9 

3.3 CONCLUSIONS 3.15 

ATTITUDE CONTROL 4.1 

4.0 INTRODUCTION 4.1 

4.1 ATTITUDE DYNAMICS BACKGROUND 4.1 

4.2 CONTROL SYSTEM REQUIREMENTS 4.2 

4.2.1 State Feedback Control Design 4.3 

4.2.2 Attitude Control Power Required 4.4 

4.2.3 Attitude Control Propellant Requirement . . .4.6 

4.2.4 Attitude Control Thruster Configuration . . 4.6 

4.3 ATTITUDE CONTROL DESIGN METHODOLOGY 4.7 

4.4 CONCLUSION 

HUMAN FACTORS 5 • 1 

5.0 INTRODUCTION 5.1 

5.1 LIFE SUPPORT SYSTEM 5.1 

5.2 ARTIFICIAL GRAVITY AND SHIP DESIGN 5.3 

5.3 RADIATION AND SHIELDING 5.6 

5.4 CONCLUSION 5.11 


V 


REPRESENTATIVE MISSION DESIGN 6.1 

6.0 INTRODUCTION 6.1 

6.1 PROCEDURE 6.1 

6.2 SAMPLE MISSIONS 6.3 

6.2.1 Saturn Envelope Mission 6.3 

6.2.2 Earth-to-Mars Mission 6.5 

6.2.3 Attitude Control System Requirements ... 6.7 

6.3 CONCLUSION 6.10 

CONCLUSION 7.2 

7 . 0 PROJECT SUMMARY 7.2 

7.1 CLOSING REMARKS 7.3 


EQUATIONS FOR AV MINIMIZATION PROBLEM A.l 


ENGINE CHARACTERISTIC CALCULATIONS B.l 


LIST OF REFERENCES R-l 


vi 


LIST OF FIGPRES 


Figure 1.1 - The Emerald City 

Figure 2.1 - The Synodic Period 2.6 

Figure 2.2 - AV Behavior Over One Synodic Period 

(Flight Time is a Parameter) 2.10 

Figure 2.3 - Assumption for Powered-Flight Phases .... 2.30 

Figure 2.4 - Geometry for Optimization Problem 2.32 

Figure 3.1 - Propellant Mass Ratio vs. 

I With AV as a Parameter 3.3 

Figure 3.2 - Impulse vs. Specific Impulse with AV as a 

Parameter and mj r . = 10 8 kg 3.4 

Figure 3.3 - Conceptual Sketch of the Gas-Core 

Nuclear Rocket 3.9 

Figure 3.4 - Number of Engines vs. AV with Powered-Flight Time 

a Parameter 3.13 

Figure 3.5 - Payload Mass Ratio vs. Engine Specific Impulse 

(I ) with AV as a Parameter 3.14 

Figure 5.1 - Radius vs. Spin Rate for Different Gravity 

Levels 5*4 

Figure 5.2 - Career Limit vs. Age 5.7 

Figure 5.3 - Method of Partial Shielding 5.8 

Figure 5.4 - Dose Rate vs. Separation for Varying Fission 

Fragment Retention Times 5.9 

Figure 6.1 - Schematic of the Emerald City 6.7 

Figure 6.2 - Non-Dimensional State Response (r = 0.5 ) . . 6.8 

Figure 6.3 - Control Torque Profile (r = 0.5 ) 6.9 

Figure 6.4 - Non-dimensional Thrust Profile (r = 0.5 ) . . 6.9 

Figure 6.5 - Sample Thruster Configuration 6.10 


vii 


LIST OF TABLES 

Table 2.1 - Eccentricity and Inclination of the Planets . . 2.8 

Table 2.2 - Percentages When Transfer is Possible Given AV 

Limit (Middle Seven Planets, Nominal Orbit = 3 AU) 2.11 
Table 2.3 - Percentages When Flight is Possible Given AV Limit 

(Inner Planets, Nominal Orbit = 3 AU) .... 2.12 

Table 2.4 - Percentages When Flight is Possible Given AV Limit 

(Inner Planets, Nominal Orbit = 4 AU ) . . . . 2.13 

Table 2.5 - Percentages When Flight is Possible Given AV Limit 

(Inner Planets, Nominal Orbit = 5 AU ) .... 2.14 

Table 2.6 - Percentages When Flight is Possible Given AV Limit 

(Inner Planets, Nominal Orbit = 10 AU) ... 2.15 

Table 2.7 - Percentages When Flight is Possible Given AV Limit 

(Outer Planets, Nominal Orbit = 3 AU) .... 2.16 

Table 2.8 - Percentages When Flight is Possible Given AV Limit 

(Outer Planets, Nominal Orbit = 4 AU) .... 2.17 

Table 2.9 - Percentages When Flight is Possible Given AV Limit 

(Outer Planets, Nominal Orbit = 5 AU) .... 2.18 

Table 2.10 - Percentages When Flight is Possible Given AV 

Limit (Outer Planets, Nominal Orbit = 10 AU) 2.19 
Table 2.11 - Synodic Periods of the Inner Planets .... 2.20 

Table 2.12 - Synodic Periods of the Outer Planets .... 2.23 

Table 2.13 - Definition of Variables Used in AV Minimization 

Problem 2.31 

Table 3.1 - Comparison of I, B , F, and m„ for LBGCNR and 


oivuv^nrt ^ v 

Table 3.2 - Engine Parameters of a SRGCNR, F^ = 4.4 x 10 5 N, 

I = 5000 s, D c = 3.3 m, and V a /V c =0.23 . . 3.11 

Table 5.1 - Possible Torus Geometries 5.4 

Table 5.2 - Continuous Dose Limits 5.7 

Table 5.3 - Solar Minimum Galactic Cosmic Ray Depth Dose 

Equivalent in Tissue as a function of Particle 

Type and LH 2 Shield Thickness 5.12 

Table 5.4 - Structural and Shield Mass Estimates .... 5.12 

Table 6.1 - Summary of Design Variables for Sample Missions 6.6 

Table 6.2 - Attitude Control System Requirements 6.8 

Table 7.1 - Envisioned Time Line for Project WISH 7.4 


viii 


NOMENCLATURE 


Greek 

|3 - thrust angle 
y - true anomaly 
A - change in 

0 - orientation from inertial reference axes 
\i - gravitational parameter, 1.3271544 x 10 11 km 3 / s 2 
x - 3.141592653 
p - density 
o - stress 

x - period, non-dimensional time 
G - orbital angular velocity 


Latin 

D - diameter 
F - force 

I - impulse, moment of inertia 
J - impulse 

N - number of crewpeople 
P - pressure, power 

R - ratio of thrust to gravity force, major radius 
TOF - time-of-f light 

V - velocity, volume 

X - inertial reference axis 

Y - inertial reference axis 
a - acceleration 

g - gravitational constant on surface of Earth, 9. 
h - height, orbital angular momentum 
m - mass 

m - mass flow rate 

n - number of engines, torus spin rate 
r - radius 
t - time, thickness 
x - state vector 


Subscripts 

RMS - root-mean-square 
T - total 
atm - atmosphere 
c - cavity 
con - for control 
dry - dry (without propellant) 
f - final state 
ff - free-flight 
i - initial , individual 
1 - payload 


of torus 


81 m/s 2 
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Subscripts (continued) 

max - maximum 
min - minor 
p - propellant 
prl - first powered flight 
pr2 - second powered flight 
sp - specific 
spe - engine specific 
syn - synodic 
th - thrust 
u - uranium 
w - engine 

0 - starting state 

1 - first state 

2 - second state 


x 





CHAPTER 1 


INTRODUCTION 


1.0 PROJECT OVERVIEW 

The year is 1620. Winter has struck with a vengeance on the 
Northeastern coast of the New World. The one hundred and two 
Pilgrims that have made the journey across the great expanse of the 
Atlantic Ocean are fighting for their lives. Because there was not 
time to establish themselves firmly on the shore, the Pilgrims are 
forced to take shelter on the Mayflower, the ship that brought them 
across the ocean, for survival. Although some will be lost, those 
that endure through the first harsh winter will do so because of 
the support of their ship. With this support, so critical during 
the initial months of the settlement, the Pilgrims will survive 
until the spring and establish a permanent colony in the unexplored 
land of North America. 

The year is now 2060. The first manned mission is arriving at 
the great ringed planet, Saturn. The purpose of this mission is to 
investigate the Saturnian moon Titan, an object of man's curiosity 
ever since the Cassini probe first explored it back at the turn of 
the century. It appears that the speculation about the possibility 
of colonizing Titan, which started as far back as the 1990's^, will 
finally be either proved conclusively or scrapped with the other 
concepts that never fulfilled their potential. Limitations of the 
interplanetary ship carrying the expedition have made cargo space 
precious. Excess amounts of food and supplies have been sacrificed 
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to accommodate the extensive scientific and data gathering equip- 
ment needed to thoroughly explore Titan and the other diverse moons 
of Saturn. Because of this, the Emerald City has taken up station 
outside of Saturn's sphere of influence to provide assistance for 
this historic mission. It will furnish food and equipment for the 
astronauts, as well as a place for them to take refuge from the 
desolate vacuum of space. It can also furnish an environment simu- 
lating that of Earth, providing artificial gravity as well as the 
psychological benefits of a large group of people, points which 
their own ship could not address. Because of the presence of the 
Emerald City, the mission can remain at Saturn for a much longer 
period of time. With the extra data that could be collected and 
recorded, future missions to the moons of Saturn may be made with 
more confidence, perhaps with the intent of establishing a perman- 
ent settlement on Titan. 

What is the connection between these two events, separated by 
over 400 years? The success of both of these scenarios is 
dependent upon the support provided by the ships in order to make 
it through the difficult initial stages of their missions. Without 
this initial support, both missions would be in serious jeopardy. 
This discussion of past and future events sets the stage for the 
introduction of Project WISH. 

1.1 PROJECT DESCRIPTION 

Project WISH (Wandering Interplanetary Space Harbor) is a 
three year design project currently being conducted at The Ohio 


1.3 


State University. This project deals with the design of a Perman- 
ently Manned Autonomous Space Oasis, designated Emerald City (see 
Figure 1.1). The purpose of Project WISH and the Emerald City is 
to serve as permanent living quarters for between 500-1000 space 
inhabitants and to provide crucial support for interplanetary 
missions in the solar system. To provide this support, the Emerald 
City will have the ability to relocate itself virtually anywhere in 
the solar system. This relocation should be accomplished in 
periods of three years or less whenever possible, thus making the 
Emerald City flexible in its mission planning. It also should have 

a lifetime of fifty years, making it an essential link in the space 

« t 

transportation network throughout the latter half of the 21 
century . 

During the 1990-1991 acedemic year. Phase II of Project WISH 
was carried out. Phase I, conducted during the 1989-1990 academic 
year, was a general level study of the major systems required for 
the Emerald City. During Phase II, a more in-depth study was con- 
ducted into the disciplines of orbital mechanics, propulsion, atti- 
tude control, and human factors. Critical elements were identified 
in all of these areas, and guidelines were established for the 
design of the Emerald City. These guidelines were then combined to 
carry out the design of two particular missions of interest, a 
Saturn Envelope mission and an Earth-to-Mars mission. The size and 
mass estimates of the Emerald City for these missions, along with 
the methods used to obtain them, are considered the main accomp- 
lishments of Phase II, and are the focus of this report. 
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CHAPTER 2 


ORBITAL MECHANICS 


2.0 INTRODUCTION 

Orbital mechanics is a central discipline in space design, and 
was one of the areas studied in-depth during Phase II of Project 
WISH. In general, orbital mechanics deals with the motion of a 
smaller body, broadly termed as the "satellite", in the gravi- 
tational field of a larger body, known as the "central body". In 
the context of Project WISH, this deals with the motion of the 
Emerald City in the gravitational field of the sun. Orbital mech- 
anics is used to find the velocity changes, or Av's, that are need- 
ed for the Emerald City to make a transfer from one place in the 
solar system to another in a certain time. Last year, some prelim- 
inary mission analysis was conducted in an ad hoc manner as an un- 
derstanding of the computer programs used in the investigation was 
developed^. A preliminary nominal orbit was chosen at 3.2 AU's, 
but this value was selected due more to the proximity of the 
asteroid belt than a detailed study of orbital mechanics. Using 
the experience gained during last year's design effort, a more 
rigorous study was performed on the transfers for the Emerald City, 
with the hope that a more substantiated nominal orbit could be 
chosen. This nominal orbit should yield the most advantageous 
transfers to as many of the planets as possible. After all, the W 
in WISH stands for "Wandering", and that is what the Emerald City 
will be doing throughout the solar system. 


2.1 NOMINAL ORBIT 


2.1.1 Background 

One of the main goals set down for Project WISH this year was 
the determination of a nominal orbit in heliocentric space. This 
nominal orbit would act as a sort of "home base" for the Emerald 
City, a place where it could return between missions for rest and 
refit. The original design requirement for the project called for 
a maximum flight time of three years from this chosen orbit to any 
planet in the solar system. The idea was to make the Emerald City 
mobile enough to be able to react to orders requiring its presence 
at a certain planet within a reasonable time frame. However, it 
was quickly discovered that this requirement placed an untenable 
burden on the prospective propulsion systems of the Emerald City, 
even extrapolating well into the 21 s ^ century. Thus the design 
goals were pared down (increasing the flight times to some planets 
while eliminating trips to other planets altogether) in order to 
make them more attainable for the time frame being considered. 

In the selection of the nominal orbit, an attempt was made to 
keep the flight times to a minimum while retaining as much of the 
solar system as possible in the "flight envelope". After much 
discussion, a new set of design goals emerged for the orbital 
mechanics and propulsion groups of Project WISH. First of all, it 
is assumed that the Emerald City will become operational sometime 
in the mid-21 st century. At this time, it is envisioned that humans 
have become firmly established on the Moon and Mars, and are begin- 
ning to look beyond. The next logical step could perhaps be human 
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exploration of the asteroid belt and the moons of Jupiter. Using 
such a scenario as background, all of the planets in the solar 
system were categorized by their perceived importance to human- 
kind's exploration of the solar system in the mid-21 s ^ century. 
These categories were employed as a guide in the decision-making 
process of what the nominal orbit should be and how much thrust the 
propulsion system of the Emerald City would need to provide. 

The first category is made up of the planets considered very 
important and consists of Earth, Mars, and Jupiter. These planets 
would more than likely be the focus of the space exploration effort 
during the next century. As mentioned earlier, colonies should 
already be established on the Moon and Mars, and intensive scien- 
tific study with preliminary exploration could be taking place on 
Jupiter and its moons. Thus, the Emerald City's ability to reach 
them in three years or less would be essential to provide support 
to these colonies and/or exploration missions. The second category 
is planets deemed important and is made up of Venus, Saturn, and 
Uranus. These are planets that would be targets for scientific 
study, but are not considered as crucial as the planets in the 
first category. Thus, it was decided that for these planets, 
transfers of up to five years would be allowable. The third and 
final category is planets considered to be of least interest, and 
includes Mercury, Neptune, and Pluto. While these planets have 
much to offer in terms of scientific value, it was felt that in the 
time frame of Project WISH, there would be little need for trips to 
these planets by a spacecraft like the Emerald City. Thus, the 
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only trips made to these planets would be in special cases when the 
planets would lie in a favorable orientation with the Emerald City. 

By using the above criteria, the design goals given at the 
beginning of last year were adjusted to more reasonable levels. 
Instead of requiring the Emerald City to make transfers from the 
nominal orbit to any planet 100% of the time in three years, this 
requirement was only applied to the planets in the very important 
category (Earth, Mars, Jupiter). For the planets in the important 
category (Venus, Saturn, Uranus), these requirements were eased 
somewhat by allowing transfers up to five years if necessary, and 
by not requiring the Emerald City to make the transfer 100% of the 
time. The design goal would instead be that the Emerald City 
should be able to make the transfer as much of the time as possible 
and still keep the propulsion system requirements somewhat 
reasonable. For the third category (Mercury, Neptune, Pluto), 
there were no specific design goals set; any transfers that would 
be possible in less than five years would be considered icing on 
the cake. Thus, with the above criteria established, the mission 
analysis was carried out. 

2.1.2 Mission Analysis 

The main difficulty in the mission analysis portion of Project 
WISH is the lack of a particular mission to analyze. By 
definition, the Emerald City is to have the capability to travel 
wherever it may be needed in the solar system within a reasonable 
time period. For example, if some crucial product is required for 


a fledgling colony on a Jovian moon, the Emerald City will not have 
the luxury of waiting years for the opening of a launch window to 
Jupiter before leaving the nominal orbit. Instead, it must be able 
to make the transfer immediately, regardless of where the planet 
might be with respect to the Emerald City. This so-called 
"transfer-on-demand" results in high Av's, which is the amount of 
velocity change that must be provided to the Emerald City in order 
for it to achieve a transfer from the nominal orbit to a planet. 
In turn, this high Av places a considerable burden on the propul- 
sion system of Project WISH. Thus, it is hoped to keep these 
mission Av's down to reasonable values in order to alleviate some 
of the propulsion problems associated with large velocity changes 
(see Chapter 3) . 

Hence, the question arises as to how to perform the mission 
analysis for the Emerald City. An ad hoc approach of picking 
target planets and launch dates for sample cases was discussed, but 
was ruled out due to its lack of rigor and tedious nature. Instead, 
a methodology was established that allowed a more thorough investi- 
gation of the transfers between a nominal orbit and the planets in 
the solar system to be performed. This methodology is based on a 
concept known as the synodic period between two orbiting bodies 
(see Figure 2.1). The synodic period is the time it takes for a 
certain angular orientation known as the phase angle between two 
bodies to repeat itself. For example, as shown in Figure 2.1, at 
some time t = 0, the Emerald City and the target planet will lie in 
a straight line connecting the two to the sun. At this point, the 


Figure 2.1 - The Synodic Period 

phase angle between the two is zero, which is known as conjunction. 
However, the Emerald City will be travelling slower than the target 
planet due to the fact that it lies further from the sun. Thus, at 
a future time, the target planet will essentially "lap" the Emerald 
City and they will once again be in conjunction at some new posi- 
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tion in their orbits. The time it takes for this to occur is the 
called the synodic period, and is denoted by t syn . The synodic 
period is related to the angular velocities of the two bodies, Qg C 
and Qj Q , and is given by 



Since the synodic period represents the time it takes for a phase 
angle to repeat itself, it can be extremely helpful in the mission 
analysis. Because the Av required to make a heliocentric transfer 
between two bodies in a given time is a function of the orientation 
of the bodies, it will be nearly cyclic in time, with a period that 
is given by the synodic period. Thus, instead of looking at trans- 
fers throughout the entire fifty year time period that the Emerald 
City is envisioned to be operating in, they can be examined over a 
single synodic period. This will result in values for AV that can 
be applied to the transfer at any future time. There are two 
assumptions on the geometry of the orbits that were made in the 
statement that the AV was cyclic with time. The first assumption 
is that the orbits of the planets are circular instead of ellip- 
tical, while the second assumption is that the nominal and plan- 
etary orbits all lie in the same plane. If these assumptions were 
invalid, then the AV would not be cyclic with time as was stated, 
but would depend on where the planet is located in its orbit. It 
can be seen from Table 2 . 1 3 that, while these assumptions are not 
perfect, they are adequate for this preliminary investigation for 
all of the planets except Mercury and Pluto. 
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Table 2.1 - Eccentricity and Inclination of the Planets 


Planet 

Eccentricity 

Incl ination 

Mercury 

0.21 

7.0° 

Venus 

0.01 

3.4° 

Earth 

0.02 

0.0° 

Mars 

0.09 

1.9° 

Jupiter 

0.05 

1.3° 

Saturn 

0.05 

2.5° 

Uranus 

0.05 

0.8° 

Neptune 

0.01 

1.8° 

Pluto 

0 . 26 

17.1° 

Eccentricity = 0 
Inclination = 0° 

==> circular orbit 
==> same orbital plane 

as Earth 


Using the computer program MULIMP* obtained from NASA Lewis 
Research Center and a simple code written to calculate the synodic 
periods for a given nominal orbit radius, the mission analysis for 
Project WISH was conducted. MULIMP is a FORTRAN code that calc- 
ulates the AV required for interplanetary transfers under the 
assumption that the change in velocity takes place instantaneously, 
i.e. infinite acceleration. Again, this assumption should not 
hamper the study because the propulsion system powered-f 1 ight times 
should be sufficiently short with respect to the overall flight 
time to keep this assumption valid for this precursory inves- 
tigation. The steps in the mission analysis then proceeded as 
follows. First, a nominal orbit radius was chosen and input into 
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the code which calculated the synodic periods for this case. Then, 
MULIMP was run for each planet at many different launch dates 
(typically between 100 and 200) over a synodic period and for 
flight times of one to five years in order to obtain the relation- 
ship between the Av, the flight time, and the fraction of time into 
the synodic period that the launch occurred at. This resulted in 
graphs similar to the one shown in Figure 2.2. This particular 
graph is for a transfer from a nominal orbit of 3 AD's (astro- 
nomical units, 1 AU is distance from Earth to Sun) to Jupiter. It 
can be seen how the Av varies over the synodic period; at some 
points it is relatively low while at other points, it is quite 
high. It can also be seen that at the completion of one synodic 
period, the AV is nearly the same as it was at the beginning, as 
expected. The small difference can be attributed to the two 
assumptions that were made about the geometries of the orbits. 
While these assumptions were made to facilitate the analysis, 
MULIMP uses the actual orbits of the planets, complete with eccen- 
tricity and inclination. Thus, the AV at each of the endpoints of 
the synodic period does not match exactly, but the fact that they 
are so close is a good indication of the validity of the assump- 
tions . 

To plot out and analyze graphs such as Figure 2.2 for all of 
the planets at each nominal orbit considered would be quite ted- 
ious. Then, the question would be what would one do with the piles 
of graphs that had been generated. To aid in the analysis of the 
data obtained from MULIMP, a program was written that would 
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Pigure 2.2 - Av Behavior Over One Synodic Period 


(Flight Time is a Parameter) 

calculate the percentage of time during the synodic period that a 
transfer from a nominal orbit to a planet would fall beneath a 
certain value of Av. Then, this data could be presented in the 
form of a chart, as shown in Table 2.2. The left-hand side of the 
chart shows what nominal orbit was used, to what planet the 
transfer was made, and what was the flight time used. For example, 
in Table 2.2, N3-P5-T3 would mean a transfer from a nominal orbit 
of 3 AU ' s to Jupiter (planet #5), with a flight time of three 
years. The numbers on the right of the chart represent the 
percentage of time during a synodic period that a transfer falls 
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Table 2.2 - Percentages When Transfer is Possible Given AV Limit 
(Middle Sev en Planets, Nominal Orbit = 3 AO) 

Case AV = 10 20 30 40 50 60 70 80 


N3-P2-T1 0.0 16.5 38.8 56.5 71.4 85.1 97.6 100 

N3-P2-T2 0.0 9.4 33.3 49.8 64.7 78.0 91.4 100 

N3-P2-T3 0.0 0.0 26.3 44.3 59.2 73.7 87.5 100 

N3-P2-T4 0.0 0.0 20.4 40.0 56.1 70.6 84.3 98.8 

N3-P2-T5 0.0 0.0 18.0 36.5 52.9 67.8 82.4 97.6 

OPTIMAL FLIGHT TIME = 1.122 

N3-P3-T1 0.0 25.6 50.2 70.0 85.0 98.2 100 100 

N3-P3-T2 0.0 28.6 48.0 64.3 78.0 91.6 100 100 

N3-P3-T3 0.0 18.1 40.5 57.3 72.7 86.8 100 100 

N3-P3-T4 0.0 8.8 34.8 52.9 68.3 83.7 99.1 100 

N3-P3-T5 0.0 0.0 31.3 49.3 65.6 81.5 96.9 100 

OPTIMAL FLIGHT TIME = 1.394 

N3-P4-T1 0.0 31.5 57.0 85.6 96.3 100 100 100 

N3-P4-T2 13.7 44.1 60.7 74.1 87.8 100 100 100 

N3-P4-T3 9.3 43.0 61.9 75.6 86.7 98.5 100 100 

N3-P4-T4 0.0 28.1 48.9 63.7 78.5 96.7 100 100 

N3-P4-T5 0.0 29.3 52.6 68.9 81.5 93.7 100 100 

OPTIMAL FLIGHT TIME = 1.764 

N3-P5-T1 0.0 7.3 20.9 29.1 34.5 40.9 92.7 100 

N3-P5-T2 11.8 32.7 83.6 99.1 100 100 100 100 

N3-P5-T3 24.5 70.0 89.1 100 100 100 100 100 

N3-P5-T4 39.1 70.0 87.3 100 100 100 100 100 

N3-P5-T5 40.9 68.2 84.5 100 100 100 100 100 

OPTIMAL FLIGHT TIME = 4.466 

N3-P6-T1 0.0 0.0 0.0 0.0 0.0 0.0 18.7 29.0 

N3-P6-T2 0.0 0.0 8.8 34.7 49.2 100 100 100 

N3-P6-T3 0.0 18.1 45.6 94.3 100 100 100 100 

N3-P6-T4 0.0 37.8 81.9 100 100 100 100 100 

N3-P6-T5 11.9 57.0 87.0 100 100 100 100 100 

OPTIMAL FLIGHT TIME = 6.860 

N3-P7-T1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

N3-P7-T2 0.0 0.0 0.0 0.0 0.0 0.0 13.8 33.5 

N3-P7-T3 0.0 0.0 0.0 0.0 27.2 48.8 100 100 

N3-P7-T4 0.0 0.0 0.0 35.4 75.2 100 100 100 

N3-P7 -T5 0.0 0.0 29.9 69.7 100 100 100 100 

OPTIMAL FLIGHT TIME = 12.960 

N3-P8-T1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

N3-P8-T2 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

N3-P8-T3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 13.7 

N3-P8-T4 0.0 0.0 0.0 0.0 0.0 18.8 42.6 91.9 

N3-P8-T5 0.0 0.0 0.0 0.0 26.9 54.3 100 100 
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below a certain Av. By using charts such as the one found in Table 
2.2, different transfer characteristics can be examined quickly and 
completely, as opposed to a piecemeal approach, which would be time 
consuming and not as thorough. Thus, now the merits of different 
nominal orbits can be compared. 


2.1.3 Analysis Results 


The resulting charts from the mission analysis can be found in 

Tables 2.3-2.10 for nominal orbits of 3, 4, 5, and 10 AO's. These 

were broken down into the inner planets in Tables 2. 3-2. 6 and the 

outer planets in Tables 2.7-2.10. 

Table 2.3 - Percentages When Plight is Possible Given AV Limit 
(Inner Planets, Nominal Orbit = 3 AO) 


Case AV = 

10 

20 

30 

40 

50 

60 

70 

80 

N3-P1-T1 

0.0 

0.0 

15.1 

31.2 

45.2 

61.3 

76.3 

90.3 

N3-P1-T2 

0.0 

0.0 

8.6 

25.8 

39.8 

54.8 

69.9 

86.0 

N3-P1-T3 

0.0 

0.0 

3.2 

21.5 

34.4 

50.5 

65.6 

80.6 

N3-P1-T4 

0.0 

0.0 

0.0 

18.3 

32.3 

46.2 

62.4 

78.5 

N3-P1-T5 

0.0 

0.0 

0.0 

15.1 

29.0 

44.1 

58.1 

76.3 

N3-P2-T1 

0.0 

16.5 

38.8 

56.5 

71.4 

85.1 

97.6 

100 

N3-P2-T2 

0.0 

9 . 4 

33 . 3 

49.8 

64.7 

78.0 

91 . 4 

100 

N3-P2-T3 

0.0 

0.0 

26.3 

44.3 

59.2 

73.7 

87.5 

100 

N3-P2-T4 

0.0 

0.0 

20 . 4 

40.0 

56.1 

70.6 

84.3 

98.8 

N3-P2-T5 

0.0 

0.0 

18.0 

36.5 

52.9 

67.8 

82.4 

97 . 6 

N3-P3-T1 

0.0 

25 . 6 

50 . 2 

70.0 

85.0 

98.2 

100 

100 

N3-P3-T2 

0.0 

28.6 

48.0 

64.3 

78.0 

91.6 

100 

100 

N3-P3-T3 

0.0 

18.1 

40 . 5 

57.3 

72.7 

86.8 

100 

100 

N3-P3-T4 

0.0 

8.8 

34.8 

52.9 

68.3 

83.7 

99.1 

100 

N3-P3-T5 

0.0 

0.0 

31.3 

49.3 

65.6 

81.5 

96.9 

100 

N3-P4-T1 

0.0 

31.5 

57.0 

85.6 

96.3 

100 

100 

100 

N3-P4-T2 

13.7 

44.1 

60.7 

74.1 

87.8 

100 

100 

100 

N3-P4-T3 

9.3 

43.0 

61.9 

75.6 

86.7 

98.5 

100 

100 

N3-P4-T4 

0.0 

28.1 

48.9 

63.7 

78.5 

96.7 

100 

100 

N3-P4-T5 

0.0 

29.3 

52.6 

68.9 

81.5 

93.7 

100 

100 
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Table 2.4 - 


Percentages When Flight is Possible Given Av Limit 
( inner Planets, Nominal Orbit 


Case AV = 

10 

20 

30 

40 

50 

60 

70 

80 

N4-P1-T1 

il 

O 11 
• II 
O II 
II 

II 

M 

O il 

O II 
II 

0.0 

tl 

II CM 
II * 
II 

II iH 
II 

33.7 

51.1 

63.0 

73.9 

N4-P1-T2 

0.0 

0.0 

0.0 

27.2 

43.5 

56.5 

66.3 

7 6.1 

N4-P1-T3 

0.0 

0.0 

0.0 

19.6 

39.1 

52.2 

64 . 1 

73 . 9 

N4-P1-T4 

o.o 

0.0 

0.0 

10.9 

3S.9 

48.9 

62.0 

71.7 

N4-P1-T5 

0.0 

0,0 

0.0 

8.7 

31.5 

46.7 

59 . 8 

69 . 6 

N4-P2-T1 

0.0 

o.o 

18.3 

41.5 

61.0 

79.3 

95.1 

100 

N4-P2-T2 

0.0 

12 . 2 

35.4 

53.7 

68.3 

82.9 

96 . 3 

100 

N4-P2-T3 

0.0 

7.3 

31.7 

51.2 

64.6 

79.3 

93 . 9 

100 

K4-P2-T4 

C .0 

0.0 

25.6 

45.1 

59.8 

74.4 

89 . 0 

100 

N4-P2-T5 

0.0 

0.0 

19.5 

41. S 

58.5 

74.4 

89.0 

100 

W4-P3-T1 

0.0 

0.0 

27 . 1 

50.7 

75.0 

93.6 

100 

100 

N4-P3-T2 

0.0 

29.3 

50.0 

66.4 

81.4 

95.7 

100 

100 

N4-P3-T3 

0 . 0 

24.3 

45.0 

61.4 

76.4 

91.4 

100 

100 

N4-P3-T4 

0.0 

17.1 

40.0 

57.1 

73.6 

87 . 9 

100 

100 

N4-P3-T5 

0.0 

10.7 

36.4 

55.0 

70.7 

85.7 

100 

100 

N4-P4-T1 

0.0 

12 . 4 

34.5 

60.2 

95.6 

100 

100 

100 

N4-P4-T2 

10 . 6 

41,6 

65 . 5 

85.8 

100 

100 

100 

100 

U4-P4-T3 

7 . 1 

40 . 7 

63.7 

77.9 

95.6 

100 

100 

100 

M4-P4-T4 

0.0 

31.9 

S4.0 

74.3 

94.7 

100 

100 

, 100 

N4-P4-T5 

0.0 

25 . 7 

54.0 

73.5 

90.3 

99.1 

100 

100 
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Table 2.5 - 


Percentages When Flight is Possibie Given AV Limit 
(Inner Planets, Nominal Orbit = 5^AU)_ 



Case AV = 

K5-P1-T1 

N5-P1-T2 

N5-P1-T3 

N5-P1-T4 

N5-P1-T5 

N5-P2-T1 

N5-P2-T2 

N5-P2-T3 

N5-P2-T4 

N5-P2-T5 

K5-P3-T1 

N5-P3-T2 

N5-P3-T3 

N5-P3-T4 

NS-P3-T5 

N5-P4-T1 

N5-P4-T2 

N5-P4-T3 

N5-P4-T4 

N5-P4-T5 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2.5 

0.0 

0.8 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

23.7 

0.0 

25 . 9 

0.0 

21.5 

0.0 

17.0 

0.0 

0.0 

0.0 

33.8 

9.4 

41.7 

0.0 

37 . 4 

0.0 

34.5 


0.0 

8.8 

7.7 

23.1 

0.0 

23.1 

0.0 

22.0 

0.0 

17.6 

0.0 

21.7 

33.3 

51.7 

31.7 

50.0 

26.7 

46.7 

24.2 

44.2 

0.0 

28.9 

46.7 

66.7 

46.7 

65.2 

43.0 

62.2 

40.0 

59.3 

12.2 

31.7 

64.0 

84.2 

64.0 

83.5 

59.7 

77.7 

56.1 

76.3 


26.4 37.4 

41.8 54.9 

40.7 56.0 

38.5 51.6 

37.4 50,5 

40.8 62.5 

68.3 84.2 

66.7 81.7 

64.2 79.2 

60 .8 77.5 

51.1 78.5 

83.7 99.3 

80.7 97.0 

77.8 94.1 

76.3 91.9 

52.5 94.2 

99.3 100 
100 100 

97.1 100 

97.1 100 


53.8 

69.2 

69.2 

83.5 

68.1 

81.3 

67.0 

79.1 

64.8 

78.0 

83.3 

100 

100 

100 

97.5 

100 

95.6 

100 

93.3 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 


Table 2.6 - Percentages When 
(Inn er Planets, 

Case Av = 10 20 30 

N10-P1-T1 0.0 0.0 0.0 
N10-P1-T2 0.0 0.0 0.0 
N10-F1-T3 0.0 0.0 5.6 
N10-Pl'T4 0.0 0.0 11.1 
N10-P1-T5 0.0 0.0 13.3 

N10-P2-T1 0.0 0.0 0.0 
N10-P2-T2 0.0 0.0 0.0 
N10-P2-T3 0.0 0.0 15.5 
N10-P2-T4 0.0 0.0 27.6 
N10-P2-T5 0.0 0.0 31.9 

N10-P3-T1 0.0 0.0 0,0 
N10-P3-T2 0.0 0.0 0.0 
N10-P3-T3 0.0 0,0 28.3 
N10-P3-T4 0.0 14.2 39.4 
N10-P3-T5 0.0 21.3 44.1 

N10'P4'T1 0.0 0.0 0.0 
N10-P4-T2 0.0 0.0 0.0 
N10-P4-T3 0.0 0.0 34,7 
N10-P4-T4 0.0 23.8 53.7 
N10-P4-T5 0.0 32.0 59.9 



40 

50 

60 

70 

I8SSSS! 

80 

0,0 

o.o 

o.o 

0.0 

14.4 

4.4 

22.2 

36.7 

51.1 

68.9 

18.9 

36.7 

53.3 

68.9 

87.8 

28.9 

43.3 

58.9 

76.7 

93.3 

30.0 

45.6 

61.1 

77.8 

94.4 

0.0 

0.0 

0.0 

o.o 

17.2 

17.2 

35.3 

56.0 

76,7 

97.4 

37 . 9 

57.8 

75.9 

96.6 

100 

48.3 

66.4 

84.5 

100 

100 

51.7 

68.1 

85.3 

100 

100 

0.0 

0.0 

0.0 

0.0 

19.7 

22.0 

43.3 

68.5 

95.3 

100 

50.4 

73.2 

95.3 

100 

100 

61.4 

81.1 

100 

100 

100 

64.6 

84.3 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

21.8 

25.9 

52.4 

87.1 

100 

100 

66,7 

94.6 

100 

100 

100 

79.6 

99.3 

100 

100 

100 

83.7 

100 

100 

100 

100 
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Table 2.7 - Percentages When Flight is Possible Given AV Limit 
(Outer Planets, Nominal Orbit = 3 AU) 


Case AV = 


N3-P5-T1 

N3-P5-T2 

N3-P5-T3 

N3-P5-T4 

N3-P5-T5 

N3-P6-T1 

N3-P6-T2 

N3-P6-T3 

N3-P6-T4 

N3-P6-T5 

N3-P7-T1 

N3-P7-T2 

N3-P7-T3 

N3-P7-T4 

N3-P7-T5 

N3-P8-T1 

N3-P8-T2 

N3-P8-T3 

N3-P8-T4 

N3-P8-T5 

N3-P9-T1 

N3-P9-T2 

N3-P9-T3 

N3-P9-T4 

N3-P9-T5 


10 

20 

0.0 

7.3 

11.8 

32.7 

24.5 

70.0 

39.1 

70.0 

40.9 

68.2 

0.0 

0.0 

0.0 

0.0 

0.0 

18.1 

0.0 

37.8 

11.9 

57.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


30 

40 

20.9 

29.1 

83.6 

99.1 

89.1 

100 

87 . 3 

100 

84.5 

100 

0.0 

0.0 

8.8 

34.7 

45.6 

94.3 

81.9 

100 

87.0 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

35.4 

29.9 

69.7 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


50 60 


34.5 

40 . 9 

100 

100 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

49.2 

100 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

27.2 

48.8 

75.2 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

18.8 

26.9 

54.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


70 

80 

92.7 

100 

100 

100 

100 

100 

100 

100 

100 

100 

18.7 

29.0 

100 

100 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

13.8 

33.5 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

13.7 

42.6 

91.9 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

19.1 

40.4 
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Table 2.8 - Percentages When Flight is Possible Given AV Limit 
(Outer Planets, Nominal Orbit = 4 AU) 


Case AV = 


N4-P5-T1 

N4-P5-T2 

N4-P5-T3 

N4-P5-T4 

N4-P5-T5 

N4-P6-T1 

N4-P6-T2 

N4-P6-T3 

N4-P6-T4 

N4-P6-T5 

N4-P7-T1 

N4-P7-T2 

N4-P7-T3 

N4-P7-T4 

N4-P7-T5 

N4-P8-T1 

N4-P8-T2 

N4-P8-T3 

N4-P8-T4 

N4-P8-T5 

N4-P9-T1 

N4-P9-T2 

N4-P9-T3 

N4-P9-T4 

N4-P9-T5 


10 

20 

0.0 

11.0 

12.3 

28.1 

22.6 

78.8 

46.6 

79.5 

56.8 

80.1 

0.0 

0.0 

0.0 

0.0 

0.0 

22.2 

2.1 

38.2 

18.8 

67.4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


30 

40 

18.5 

24.0 

50.7 

99.3 

95.9 

100 

94.5 

100 

93.8 

100 

0.0 

0.0' 

20.1 

32.6 

45.1 

99.3 

94.4 

99.3 

96.5 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5.4 

36.5 

31.1 

73.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


50 

60 

29.5 

36.3 

100 

100 

100 

100 

100 

100 

100 

100 

1.4 

16.7 

41.7 

98.6 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

28.4 

42.6 

70.9 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

20.2 

27.9 

49.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


70 80 


43.8 

94.5 

100 

100 

100 

100 

100 

100 

100 

100 

20.8 

25.0 

100 

100 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

20.3 

29.7 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

17.1 

41.1 

94.6 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

15.4 

41.0 


2.17 


Table 2.9 - Percentages When Flight is Possible Given AV Limit 
(Outer Planets, Nominal Orbit = 5 AU) 


Case AV = 

N5-P5-T1 

N5-P5-T2 

N5-P5-T3 

N5-P5-T4 

N5-P5-T5 

N5-P6-T1 

N5-P6-T2 

N5-P6-T3 

N5-P6-T4 

N5-P6-T5 

N5-P7-T1 

N5-P7-T2 

N5-P7-T3 

N5-P7-T4 

N5-P7-T5 

N5-P8-T1 

N5-P8-T2 

N5-P8-T3 

N5-P8-T4 

N5-P8-T5 

N5-P9-T1 

N5-P9-T2 

N5-P9-T3 

N5-P9-T4 

N5-P9-T5 


10 20 


6.2 

15.8 

15.8 

26.6 

20 . 9 

57.1 

32.8 

87 . 6 

64.4 

84.7 

0.0 

0.0 

0.0 

4.8 

0.0 

23.8 

11.6 

36.7 

20.4 

73.5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


30 

40 

18.1 

22 . 6 

40.7 

100 

98.9 

100 

99.4 

100 

98.9 

100 

0.0 

0.0 

22.4 

29.3 

40.1 

98 . 6 

100 

100 

99.3 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

8.2 

35.4 

31.6 

65.2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


50 60 


24.3 

29.4 

100 

100 

100 

100 

100 

100 

100 

100 

12.9 

16.3 

36.1 

97.3 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

27.8 

38.0 

54.4 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

23.1 

29.3 

47.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


70 80 


35.0 

40.7 

100 

100 

100 

100 

100 

100 

100 

100 

19.7 

22.4 

100 

100 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

20.3 

27.2 

100 

100 

100 

100 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

21.1 

38.1 

60.5 

100 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

18.0 

39.5 
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Table 2.10 - Percentages When Flight is Possible Given AV Limit 
(Outer Planets, Nominal Orbit = 10 AU) 


Case AV = 

10 

20 

30 

40 

50 

60 

70 

80 

N10-P5-T1 

0.0 

0.0 

0.0 

0.0 

11.0 

14.8 

17.4 

20.6 

N10-P5-T2 

0.0 

0.0 

20.6 

27.1 

32.9 

41.9 

100 

100 

N10-P5-T3 

0.0 

21.3 

36.1 

67.1 

100 

100 

100 

100 

N10-P5-T4 

7.7 

35.5 

94.8 

100 

100 

100 

100 

100 

N10-P5-T5 

17 . 4 

54.2 

100 

100 

100 

100 

100 

100 

N10-P6-T1 

5.2 

7.0 

7.4 

10.0 . 

13.0 

13.9 

14.3 

17.8 

N10-P6-T2 

6.5 

11.3 

14.8 

19.1 

22.2 

27.8 

33.0 

41.7 

N10-P6-T3 

8.3 

20 . 4 

27.0 

33.5 

44.8 

100 

100 

100 

N10-P6-T4 

15.7 

27.4 

35.2 

97.0 

100 

100 

100 

100 

N10-P6-T5 

19.1 

34.3 

97.0 

100 

100 

100 

100 

100 

N10-P7-T1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

N10-P7-T2 

0.0 

0.0 

0.0 

0.0 

7.8 

12.4 

15.7 

19.0 

N10-P7-T3 

0.0 

0.0 

0.0 

15.0 

19.6 

24.8 

30.1 

36.6 

N10-P7-T4 

0.0 

0.0 

16.3 

22.9 

29.4 

38.6 

100 

100 

N10-P7-T5 

0.0 

10.5 

23.5 

32.0 

87.6 

100 

100 

100 

N10-P8-T1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

N10-P8-T2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

N10-P8-T3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

15.5 

20.8 

N10-P8-T4 

0.0 

0.0 

0.0 

0.0 

15.5 

22.6 

29.2 

36.3 

N10-P8-T5 

0.0 

0.0 

0.0 

16.7 

26.2 

34.5 

47.6 

100 

N10-P9-T1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

N10-P9-T2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

N10-P9-T3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

N10-P9-T4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 

N10-P9-T5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

20 . 9 

29.9 
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Inner Planets 


After reading the above discussion and examining Tables 2.3- 
2.6, one might be tempted to arrive at a conclusion similar to the 
following. "Well, since planet Earth fell into the first category 
of planets that had to be reached 100% of the time in three years 
or less, and according to Tables 2. 3-2. 6, Earth can be reached 100% 
of the time in three years for Av's greater than 60 km/ s , then the 
propulsion system will have to provide a one-way Av of greater than 
60 km/s for the design goals to be met." However, upon closer 
inspection, this conclusion is incorrect due mainly to two reasons. 
The first reason is that the synodic period of the inner planets 


are rather low. 
Table 2.11 - 

as can be seen 
Synodic Periods 

in Table 2.7. 
of the Inner 

They range 
Planets 

from about 

Planet 

Mercury 

Venus 

Earth 

Mars 

T srn (years) 


- 



n = 3 AU 

. 25 

.70 

1.24 

2.95 

n = 4 AU 

. 25 

. 67 

1.14 

2 . 46 

n = 5 AU 

. 25 

. 65 

1.10 

2.26 


one-quarter of a year for Mercury to over two years for Mars. This 
is a key point in determining what Av's will need to be provided 
for transfers to these inner planets. Because the synodic period 
is so short, the Emerald City would pass through the entire period 
in less than three years time. Thus, the Emerald City would have 
the luxury of waiting at its nominal orbit for some amount of time 
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after receiving its orders to move, until the Av requirements are 
sufficiently low to make a transfer possible in less than three 
years and it would still arrive within three years of the time it 
received its orders to go somewhere. 

For example, say that the Emerald City is in a nominal orbit 
at 4 AU's from the sun. On February 10, 2051, it receives orders 
to travel to Earth, and its arrival is required within two years. 
On February 10, the Av required to make a two year transfer at this 
time is found to be 63.9 km/s, a rather high value that would place 
a large burden on the propulsion system. So instead, the Emerald 
City decides to wait at the nominal orbit until July 3, 2051, when 
the AV required to make a one year transfer is only 23.1 km/sec. 
This value is much more easily obtained by the propulsion system, 
and the Emerald City arrives at Earth on July 2, 2052, about 222 
days ahead of the time it was required to be there. This scenario, 
calculated using MULIMP, is possible because the synodic period for 
Earth is only 1.14 years, short enough to allow all of the possible 
phase angles to present themselves to the Emerald City within the 
time that it is required for it to make the transfer. Thus, 
instead of leaving immediately, it can wait for some length of 
time, and then make a transfer that, even with a shorter flight 
time, has a lower AV because of a more favorable orientation. 

The second reason that the percentages that a transfer is 
possible for a given flight time cannot be taken at face value is 
because of the following fact. In this analysis, the flight time 
for each transfer is specified, and the AV required to make this 
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exact transfer is found. One might expect that this AV would go 
down as the flight time is made longer, because of the less 
demanding duration. In general, however, this is not the case for 
the inner planets. The optimal flight time (the flight time 
resulting in the lowest AV) for the inner planets is on the order 
of one to two years. Thus, by forcing a trip to one of these 
planets in, say, three years, the AV required actually increases. 
This is because the trajectory will not go directly to the target 
planet, but will stall by taking an indirect route in order for the 
flight time to be the specified value. This phenomenon can be seen 
in the following example, predicted by MULIMP. The Emerald City 
receives orders to proceed to Venus on February 27, 2050, and is 
ordered to be there within 2.5 years. At this point, if the 
Emerald City were to try to make a transfer with a flight time of 
2.5 years, the AV required would be 46.1 km/s. However, by making 
the transfer in only 2.33 years, the AV would only be 18.8 km/s, a 
much better result. 

Thus, because of the nature of the transfers to the inner 
planets for the two reasons cited above, the actual percentages of 
time that a flight is possible for a given AV is rendered 
meaningless. Instead, what should be looked at is the lowest AV a 
transfer can be made at all with a flight time of three years minus 
the synodic period. If a transfer can be made at any point in the 
synodic period for this flight time, then the transfer can be made 
100% of the time. Again using the Earth example, in 1.14 years, 
all of the phase angles between it and the Emerald City will have 
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occurred due to the length of the synodic period. Thus, if Earth 
can be reached any percent of the time in 1.86 years or less with 
a given AV, then it can be said that the transfer can be made 100% 
of the time in under three years from receiving its orders, because 
it can wait until the favorable orientation appears. 

♦ 

Outer Planets 

The results from the mission analysis for the outer planets 
(Jupiter to Pluto) can be found in Tables 2.7-2.10 for nominal 
orbits of 3, 4, 5, and 10 AU's. The investigation of the transfers 
to the outer planets must be treated differently than to the inner 
planets because the synodic periods are much longer for the outer 
planets, as seen in Table 2.12. Because of this, the Emerald City 
Table 2.12 - Synodic Periods of the Outer Planets 


Planet 

Jupiter 

Saturn 

Uranus 

Neptune 

Pluto 

T syn (years) 






n = 3 AU 

9.26 

6.31 

5.54 

5 . 37 

5.31 

n = 4 AU 

24.62 

10.98 

8.85 

8.41 

8.27 

n = 5 AU 

196.66 

18.02 

12.91 

12.00 

11.71 


will encounter only a fraction of the possible phase angles, and 
must be prepared to make transfers at less than optimum conditions. 
Thus, for the outer planets, the percentage of time that a transfer 
is possible for a fixed AV becomes important, a yardstick that 
different nominal orbits can be measured by. 
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It should be noted what effect moving the nominal orbit 
further from the sun had on the ability of the Emerald City to 
reach the outer planets. As seen in Tables 2.7-2.10, there was 
little benefit gained by moving the nominal orbit away from the 
sun, and in some cases, the percentages actually went down. This 
was an important consideration for the selection of the nominal 
orbit . 

2.1.4 Nominal Orbit Selection 

The sample nominal orbits that were examined were at 3, 4, 5, 
and 10 AO's. After examining the percentage tables generated by 
MULIMP and some discussion, an orbit at 4 AU's was selected as the 
preliminary nominal orbit. This selection was based on several 
factors. The first one is the central location of this orbit 
between Mars, at 1 . 5 AU and Jupiter, at 5.2 AU. This orbit puts 
the Emerald City relatively near to the asteroid belt, where vast 
deposits of valuable materials can be found, and close to the 
presumed focal points of interplanetary space activity in the mid- 
21 s ^ century. Mars and Jupiter. From this orbit, a three-year 
transfer to Jupiter can be accomplished 100% of the time with a 
one-way AV of 40 km/s. Thus, if the goals that were laid down at 
the beginning of this chapter are to be met, the propulsion system 
would have to provide at least this much AV capability. Also, 
Venus, Earth and Mars can be reached under AV's of 40 km/s in three 
years or less, considering the points that were brought up in the 
discussion of the inner planets previously. Saturn can be reached 
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nearly 100% of the time at 40 km/s using flight times of three 

years, which is an added bonus. To reach Uranus with a Av of 40 

km/s, flight times of four and five years must be used, and, for 

five years, the transfer is possible only 73% of the time. This 

was felt to be acceptable, because Uranus fell into the category of 

s t 

planets deemed marginally important for exploration in the mid-21 
century. Mercury can be reached 100% of the time if a trip to that 
planet is required, but Neptune and Pluto cannot be reached at all 
by transfers with flight times up to five years and a Av limit of 
40 km/s. Thus, trips to these outer planets would require special 
considerations, such as somehow refueling part of the way there or 
increasing the propulsion capabilities. A nominal orbit of 4 AU's 
also puts the Emerald City within reasonable distance of Earth for 
communications and emergencies, which is why a larger radius orbit 
was not selected (the benefits gained by moving the orbit out were 
generally negligible). 

An interesting point that should be brought up when discussing 
the necessary Av that the propulsion system must provide is the 
return to the nominal orbit. Early in this study, it was generally 
assumed that accounting for the return trip would roughly double 
the maximum AV that would be necessary. However, upon further 
investigation, some doubts were raised as to the validity of this 
assumption. When using the MULIMP code, the orbit of the Emerald 
City is entered by using the six classical orbital elements; the 
eccentricity, the semi -major axis, the inclination, the longitude 
of the ascending node, the argument of the periapsis, and the time 
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of periapsis passage. For the circular orbits that were examined, 
there is no periapsis, and the last element is replaced by some 
value that locates the Emerald City in the orbit at a given time. 
When going from the nominal orbit to a target planet, the Emerald 
City will leave from a particular point on the nominal orbit and 
arrive at the target planet, and the AV can be calculated by 
MULIMP. However, for the return trip, there is nothing that says 
the Emerald City has to return to a particular point on the nominal 
orbit. Instead, it is free to return to any point along the 
nominal orbit. Unfortunately, when the orbital elements are used 
to define the nominal orbit, a "ghost point" is created that moves 
around the defined nominal orbit, as if it were a planet in that 
orbit. If MULIMP were used to figure out a AV required for the 
return trip, it would calculate it as if it had to return to the 
"ghost point's" location instead of anywhere on that orbit. 
Because of this freedom, the Emerald City could return to the 
nominal orbit at the point that would require the least AV for the 
transfer, and the AV's would probably not be as high as originally 
expected. It was not known if or how MULIMP could be used to 
calculate the AV for a situation like this, so an total AV of 50 
km/ s was assumed as the necessary upper limit to give to the 
propulsion system. For the inner planets, this figure should be 
quite adequate due to the Emerald City's ability to make optimal 
transfers as was discussed above. For the outer planets, it was 
felt that this figure should result in 100% ability to reach 
Jupiter and Saturn in three years and above 50% to reach Uranus in 
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five years (or less, in some cases). Of course, Neptune and Pluto 
still are unreachable without some special considerations. 

2.2 AV MINIMIZATION 
2.2.1 Introduction 

The field of optimization is a powerful tool in the 
engineering design process. The basic optimization problem is one 
of determining the minimum or maximum a function, called the 
objective function, can obtain, subject to certain constraints on 
its variables. This problem is represented by 

min^gjj* f(x) (2.2) 

subject to the constraints 

g^(x) = 0 for j = 1 ,...., m e 
g^(x) > 0 for j = m e +l,....,m 

Xj < X < X u . 

The function f(x) is the objective function and x is the design 
variable vector of length n, i.e. n design variables. There are m 
constraints, the first type being called equality constraints while 
the second type is called inequality constraints. The third type 
of constraint is on the design variables themselves and specifies 
the range of values that these variables can take. There has been 
much theory developed on techniques to solve optimization problems 
in the above form. The particular method selected depends upon the 
nature of the problem; the type of constraints, whether the problem 
is linear or non-linear, whether the gradients of the functions are 
to be eval uated anal yti cal 1 y or by finite-difference methods , etc.. 
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All of these factors will need to be considered when selecting the 
computer code to be used to solve the optimization problem. 

With the above discussion in mind, an optimization problem 
studied this year will be discussed, namely the problem of Av 
minimization . 

2.2.2 Background 

One of the major problems that has surfaced in the first two 
years of development of Project WISH deals with the transfers from 
a nominal orbit to the planets in the solar system. It has been 
noted that the Av needed for these transfers have a great impact on 
the propulsion system requirements. As the values for the AV 
increase, the demands on the propulsion system become considerable. 
Besides straining the propulsion system, this has a cascade effect 
on the rest of the systems of the Emerald City; greater propulsion 
system requirements (thrusts, powered-flight times...) lead to more 
difficult problems in shielding, structural dynamics, propellant 
storage, etc... Thus, it is in the best, interest of the Emerald 
City to keep these AV * s as low as possible. For this reason, it 
was felt that an optimization of Av was justified in order to keep 
propulsion system requirements as reasonable as possible. It may 
seem that this would be redundant with the work that was presented 
earlier in this chapter. However, the analysis presented there 
used a computer program that assumed an instantaneous change in 
velocity, with no regard for the powered-flight times, thrust 
required, or the thrust angles that could be used. By performing 


2.28 


this optimization, it is hoped that these factors can be taken into 
account, as well as determining the orientation between the Emerald 
City and the target orbit, typically near the destination planet, 
that will yield the lowest AV. 

2.2.3 Assumptions 

In order to carry out the optimization, certain assumptions 
regarding the nature of the orbits and transfers had to be made. 
The first assumption is that planets are all in circular orbits and 
are all in the ecliptic plane, i.e. the same plane as the Earth. 
As was shown in Table 2.1, this is a reasonable assumption for all 
of the planets except Mercury and Pluto, which are of the least 
interest in Project WISH. The second assumption made was that 
there is no loss of propellant during the free-f light phase of the 
transfer. Of course, there will be some boil-off of the liquid 
hydrogen due to heat leakage through the storage tanks, but for 
simplicity, it was ignored for this analysis. The validity of this 
assumption is questionable, but it was reasoned that the results 
would still be applicable for this preliminary study. By making 
this assumption, the acceleration at the end of the first powered- 
flight phase was set equal to the acceleration at the beginning of 
the second powered-flight phase. 

The third, and most pivotal, assumption is that the trajectory 
during the powered-flight times is unaffected by the pull of the 
sun’s gravity (see Figure 2.3). During the powered-flight times, 
the Emerald City is assumed to follow a pro jecti 1 e- 1 ike course. 
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Figure 2.3 - Assumption for Powered-Flight Phases 


determined by the X- and Y-components of the initial velocity and 
acceleration (thrust) vectors. This assumption is based on the 
value of the ratio of accelerations due to thrust and gravity 
during the powered-flight phases, given by 


m{t) 



\un(C) 


- m -tilt n -1.327 X 10 11 

s 2 


(2.3) 


It can be seen that R is dependent upon the thrust, , the radius, 
r(t), and ship mass, m(t). For representative values for the 
Emerald City, R will be a very large number, and thus it was felt 
that this assumption was valid for this preliminary investigation. 
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2.2.4 Geometry 


The geometry that was used for the Av optimization is shown in 
Figure 2.4. The technique used to model the interplanetary 
transfer was to break down the transfer into three arcs. The first 
and last arcs represents the initial and terminal powered-flight 
phases, respectively, while the middle arc represents the free- 
flight phase. In the powered-flight portions, the thrust level, 
powered-flight times, and thrust angles govern the trajectory, 
while in the free-flight portion, the classical relationships 
governing orbiting bodies are used (more on these relationships can 
be found in Appendix A). The variables shown in Figure 2.4 are 
defined in Table 2.13. 

Table 2.13 - Definition of Variables Used in 
AV Minimization Problem 

r fl - nominal orbit radius 

r^ - radius at end of first powered 

flight, beginning of free-flight 

r 2 - radius at end of free-flight, 
beginning of second powered 
flight 

r^ - target orbit radius 

9g - angle between inertial reference frame 
X-axis and Tq vector 

9^ - angle between inertial reference frame 
X-axis and r^ vector 

“ angle between inertial reference frame 
X-axis and r? vector 

0^ - angle between inertial reference frame 
X-axis and r^ vector 

- angle between inertial reference frame 
X-axis and target planet at launch 

Yj - angle between inertial reference frame 
X-axis and eccentricity vector ( line 
of apsides) 

y, - true anomaly of free-flight orbit at 
point 1 

Yj ~ true anomaly of free-flight orbit at 
point 2 
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2.2.5 Methodology 


In an optimization problem, there is a objective function that 
represents the variable to be optimized and constraint functions 
that represent the restrictions on the problem. For the AV 
minimization problem, the general objective function is 

AV r - AVi + AV^ - (Vi - V a ) + ( V f - V 2 ) (2.4) 

where V^, V Q , Vj , and Vj are the velocities at points 1, 0, f, and 
2, respectively (see Figure 2.4). The general constraints for this 
problem are given by constraints on the total time-of-f light, the 
acceleration levels, the powered-flight time lengths, the 
rendezvous requirements, and physical variable constraints. The 
acceleration and powered-flight-time constraints arise from the 
limitations of the propulsion system, while the total time-of- 
f light constraint governs the allowable time to make the transfer. 
The rendezvous constraint requires that the position and velocity 
of the Emerald City at the end of the transfer matches the position 
and velocity of the target orbit near the destination planet at 
that time, and the final constraint is to ensure that there are no 
physically meaningless solutions found, such as negative times or 
negative magnitudes. Thus the problem is to find the equations 
that relate Vg , v^ , Vj , and v^ to the design variables chosen and 
determine equations for the constraints listed above. 

The first step that was taken was to define a state vector for 
point i, x^, such that 

*i (2.5) 
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The values at the end of each powered-flight phase can be related 
to the values at the beginning of each powered-flight phase by 


- X 0 + 

X f “ *2 + 


( 2 . 6 ) 


where 


LXj-iAfj L (2.7) 

By examining the equations governing the trajectory of the Emerald 
City during the powered and free-f light phases, and by contem- 
plating the nature of the problem, eight design variables were 
identified. These variables are the two powered-flight times, t ?r j 
and tp r2 , the two thrust angles, and |3 2 , acceleration level, 
a 0 , the true anomaly at the end of free-flight, y 2 , and the initial 
orientations of the Emerald City and the target planet, and 6^, 
respectively. The functional dependencies between these variables 
and the state vectors were found to be 


x 1 - X 0 (d 0 ) + A*! ^ 

Xf - X 2 {r 2 )\g i + Aj2j(y 2 ,a o' tprl' ^pr2^ 

where the vertical bar signifies dependance on a previous state. 
After applying equations for and AV 2 as found in Appendix A, the 
objective function was written as 

AV T - -J_g(ln(l - iS^EEi) + ln(l - a o t P r 2 ) } (2 .9) 

^ i ep sr I,*# - a otpTi 
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where I jp is the specific impulse of the propulsion system (see 
Chapter 3), and g is the gravity constant on the surface of the 
Earth, 9.81 m/s^. The constraint functions may then be written as 

1- tff + tpri + tp r 2 < T0F mi 

2 - a o < a o,.ax 

3 • fc prl + fc pr2 < Sr-mai 

4 ‘ x f,EC = x f , TO * 

a o' t prl ' fc pr2' > 

It turns out that the objective function is dependant only on three 
of the eight design variables, namely a fl , t pr j, and t pr j . The other 
five design variables enter into the problem through the first and 
fourth constraints above. These constraints on the total time-of- 
flight and the rendezvous requirement incorporate the rest of the 
design variables through the equations governing the powered and 
free-flight phases of the trajectory. These equations can be found 
in detail in Appendix A. 

Once the objective function and constraint functions are 
determined, an optimization code can be utilized to solve this 
problem. A potential optimization subroutine was identified, but 
due to time and personnel constraints, the problem was not 
programmed. This subroutine solves a general, non-linear 
optimization problem using a successive quadratic programming 
algorithm and finite difference gradients. It is hoped that this 
optimization problem can be programmed during next year's design 
effort and some meaningful results obtained. 
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2.3 CONCLUSIONS 


This year, a detailed study was performed on the orbital 
mechanics of Project WISH. Based on this study, a nominal orbit of 
4 AU’s was selected, which places the Emerald City roughly 
600,000,000 kilometers from the sun, between the orbits of Mars and 
Jupiter. From this orbit, it was determined that the planets 
Mercury, Venus, Earth, Mars, Jupiter, and Saturn can be reached 
100% of the time in flight times of three years or less with a Av 
limit of 40 km/ s one-way. Uranus can be reached 73% of the time 
with this AV in flight times of five years, which was felt to be 
adequate for this project. Neptune and Pluto cannot be reached at 
all at this AV with flight times of up to five years, but it was 
believed that there would be little need for a ship such as the 
Emerald City to make trips to these planets during the time frame 
considered. The AV minimization problem was also examined and 
equations were generated that yield the trajectory of the Emerald 
City during both the powered- and free-f light phases. The 
objective and constraint functions were determined, but due to a 
lack of time, the problem was not programmed into an optimization 
code. It is hoped that in the future, this problem can be solved 
and meaningful results can be obtained. Also, more work should be 
done on determining the Av's required to return to the nominal 
orbit from a planet. With additional study in this area, a more 
accurate round-trip figure for the Av's can hopefully be obtained. 
As of the end of this study, a total round-trip AV of 50 km/ s was 
used as the upper ceiling for the Emerald City. 
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CHAPTER 3 


PROPULSION 


3.0 INTRODUCTION 

The propulsion system of the Emerald City plays an important 
role in the success of Project WISH. For this reason, it was one 
of the most rigorously studied subsystems this year. At the 
conclusion of the Phase I design conducted last year, an antimatter 
engine was selected as the main propulsion system of the Emerald 
City^. However, it was decided during Phase II that antimatter was 
too theoretical for the time frame of Project WISH, and that the 
propulsion system chosen should be more readily procurable by the 
mid-21 s ^ century. With the hope of bringing the choice of the 
propulsion system into the Emerald City's time frame, the 
propulsion team looked at several areas of interest in spacecraft 
propulsion. First, a general level study was conducted in order to 
determine some of the important parameters associated with 
spacecraft propulsion and what impact they had on the choice of a 
particular system. Once this was accomplished, a propulsion system 
was chosen and a specific propulsion study was conducted in order 
to determine the significant characteristics of this system. From 
this study, it was possible to calculate other important parameters 
that were dependent upon the choice of a particular propulsion 
system. Finally, this propulsion system was used in the design of 
two representative missions, which may be found in Chapter 6. 
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3.1 GENERAL PROPULSION STUDY 


3.1.1 Theoretical Background 

The purpose of this general level study is to determine 
which parameters play a major role in the choice of a propulsion 
system, and what the orders of magnitude of these parameters need 
to be for an application such as Project WISH. The study is 
independent of the type of engine that will be used for Project 
WISH; only the physical nature of spacecraft propulsion enters into 
this investigation. The first step in the analysis was to deter- 
mine the effects that the mission requirements have on propulsion 
parameters. Since the Emerald City is expected to be a large 
spacecraft, the propulsion system will be required to produce high 
impulse levels, J, at high specific impulses, I sp . High specific 
impulse, which is the impulse produced per unit weight of 
propellant expelled, is desirable because it corresponds to less 
propellant required to accomplish a mission. The I sp plays an 
important role in the equation for the propellant mass ratio, given 
by 

- i - (3.1) 

which is plotted in Figure 3.1. This figure demonstrates how the 
propellant mass ratio, which is the fraction of the total initial 
ship mass consisting of propellant, decreases as the I sp is 
increased or the AV is decreased. Thus, to keep the propellant 
mass required for a given mission down, a propulsion system capable 
of generating a high I jp should be utilized. 
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However, for the Emerald City, a high-Ijp engine alone is not 
enough. An engine capable of producing high impulse levels as well 
as the high I sp 's is needed. The impulse required for a mission is 
a function of the dry mass of the ship (no propellant included) and 
the A V and is given by 

-All r> 

J- m^I^g {e 1 -* - 1) (3 ' 2) 

which is plotted in Figure 3.2. From this figure, it can be seen 
that the impulse required for an Emerald City mission will be very 
large, on the order of 10 13 to 10 14 newtons. This is due to two 
factors; the large dry mass of the Emerald City, and the high AV 
levels that will be required for its missions (see Chapter 2). The 
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effect that the I. has on the impulse required can also be seen 

sp 

from Figure 3.2. As the I, n goes down, the impulse required goes 
up, eventually approaching infinity as I sp approaches zero. This 
is because of the extra propellant mass that must be carried at the 
lower I 's, which makes the Emerald City heavier and more difficult 

sp 

to move. Since impulse is the product of a force and a time, high 

impulses would necessitate high thrusts, F, and high powered-flight 

times, t... Neither of these are desirable from an engine 
pr 

maintenance/ 1 i fetime viewpoint. It may also be seen from Figure 
3.2 that, if the Emerald City were to burn for its entire flight 
time (~ 3 years or ~ 10® seconds ), the thrust that would be 
required would still be on the order of 10® newtons. 
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Once it has been determined that high Ijp's and high thrusts 
are needed, the next step was to search for a propulsion system 
that could satisfy these requirements. A general level study was 
conducted in an attempt to find such a system from choices that 
included chemical, electrical, nuclear, and antimatter. However, 
antimatter had already been ruled out for being too theoretical in 
nature, even for an advanced project such as Project WISH. Thus, 
the remaining systems were analyzed in terms of the specific 
impulse and thrust levels that each system could provided. 

3.1.2 General Propulsion Systems 

Chemical propulsion was the first type of propulsion analyzed 
and was promptly rejected. Although the thrust values were 
adequate, the specific impulses of these systems fell far short of 
the requirements that the Emerald City would place on them. Elec- 
trical propulsion was the next system that was studied. These 
systems produced high enough Ijp's, but were severely limited in the 
amount of thrust that could be produced. Nuclear thermal pro- 
pulsion ( NTP ) was the final type of system examined. Estimates of 
the nuclear propulsion systems of the 21 s ^ century gave the best 
combination of high thrust and high I jp that would be required for 
Project WISH. Thus, NTP was examined in more detail as a primary 
propulsion system for the Emerald City. 

One of the main points considered when examining nuclear 
systems was the question of fission or fusion. The propulsion 
group decided that fusion, like antimatter, would be too 
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theoretical for the Emerald City. Another drawback of using fusion 

propulsion is the mass penalty that would accompany the system. 

Compared to the fission systems, the fusion systems are estimated 

to be substantially more massive. Of course, a breakthrough in 

fusion research could lead to the development of an operational 

fusion system suitable for the Emerald City, but for now, fusion 

was rejected as a main propulsion candidate. Fission systems have 

the advantage of being more advanced and currently less massive 

than their fusion counterparts. Unfortunately, associated with 

fission is a rather troublesome radiation problem, which leads to 

an increase in the shielding mass required for human safety. Also, 

the more developed types of fission systems have relatively low 

I *s for this project, and the more advanced engines witli the 
sp 

higher I 's are still conceptual in nature. With all of these 

sp 

considerations in mind, it was decided that a fission nuclear 
propulsion system would be used for the Phase II design of the 
Emerald City. 

Within the fission category of NTP, there are two major 
subgroups; solid-core reactors and gas-core reactors. As its name 
suggests, a solid-core NTP system utilizes a fissioning solid-core 
of uranium to generate thermal energy. This energy is transferred 
to a lightweight propellant, usually hydrogen, which is then 
expelled through a convergent-divergent nozzle. However, the core 
is limited in temperature by material considerations, and thus the 
achievable I sp is similarly limited. Current estimates and some 
testing during the NERVA (Nuclear Energy Rocket Vehicle Appli- 
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cation) place the solid-core I.„ in the range of 800-1100 seconds*. 

•>r 

The gas-core NTP system attempts to overcome this issue of limited 
temperature by allowing the core temperature to increase beyond the 
boiling point of uranium, thus resulting in a gaseous core. This 
increase in temperature enables more thermal energy to be trans- 
mitted to the propellant, which results in a higher possible I jp . 
Current theoretical studies give a possible I Jp for this conceptual 
engine on the order of 2000-8000 seconds*. It was felt that the 
I $ p of the solid-core rockets were simply too low for the needs of 
Project WISH (see Figures 3.1 and 3.2) and thus the gas-core NTP 
system, with its higher estimated attainable thrusts and I sp 's was 
chosen for further study as the main propulsion system of the 
Emerald City. 

There were two major types of gas-core nuclear rockets (GCNR) 
studied, the light bulb (closed-cycle, LBGCNR) and space radiator 
(open-cycle, SRGCNR) . The closed-cycle gas-core engine has its 
gaseous core encased in a transparent material that allows radia- 
tion from the fissioning core to heat the propellant while keeping 
the propellant and core physically apart. The open-cycle gas-core 
engine has the propellant and core in direct contact, which allows 
for better heat transfer, but also leads to a loss of uranium from 
the core. In determining which of these propulsion systems was 
best for the Emerald City, three main engine parameters were 
investigated, the I sp , the thrust per engine, , and the mass per 
engine, . Table 3.1 shows some comparison figures of these three 
parameters for the two systems considered*. This table shows that 
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the closed-cycle engine has a significantly lower specific impulse 
than the open-cycle engine. However, the maximum expected thrust 
value for the closed-cycle system is higher than that of the open- 
cycle engine . It will Table 3.1 - Comparison 


be shown later in this 
chapter in Figure 3.5 
that the advantage of 
producing higher thrust 
decreases as the thrust 
level increases. It was 
felt that the closed- 
cycle gas-core engines 
still had I,„ values that 
were too low for use by 
the Emerald City, and 


LBGCNR DATA 


Thrust (N) 

m^ (kg) 

(sec) 

133,370 

14,050 

1780 

1,334,200 

34,475 

2355 

4,002,800 

385,500 

2635 1 


SRGCNR DATi 


Thrust (N) 

m^ (kg) 

Icn (sec) 

22,240 

36,280 

2400 

177,900 

101,440 

5500 

444,750 

213,350 

6000 


thus, the open-cycle engine (hereafter referred to as SRGCNR) was 
chosen as the propulsion system for the Emerald City. The values 
of thrust and I S(J that were used as representative figures for the 

7 

SRGCNR in the remainder of this report are 1 

I $ p = 5000 seconds 
F i = 4.4 x 10 S N. 

These values fell within the range of thrusts and I sp ' s that were 
discussed in Reference 7, and were felt to be reasonably obtainable 
if the SRCGNR is ever developed. Once the propulsion system and 
important parameters were determined, it was then possible to con- 
duct a detailed study of the SRCGNR system of the Emerald City. 
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3.2 SPACE RADIATED GAS-CORE NUCLEAR ROCKET STUDY 

With the specific propulsion system determined and values for 
I S p and F^ known, it was then possible to perform a more detailed 
study of this system. Because the SRCGNR is so conceptual in 
nature, there are many unanswered questions about its form and 
function. Figure 3.3 is a schematic of the GCNR that was con- 
ceptually developed at NASA Lewis Research Center by Mr. Robert 
Ragsdale^. The main mechanical components of this engine are the 
pressure shell, the moderator, the turbopump, the nozzle, and the 


Radiator 



Figure 3.3 - Conceptual Sketch of the Gas-Core Nuclear Rocket 
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radiator, which work together in the following manner. First, 
liquid hydrogen is pumped through the moderator by the turbopump to 
aid in heat removal there. Then, the hydrogen is pumped into the 
cavity where the uranium isotope, U^, is undergoing the fission 
process. The hydrogen absorbs heat from the fissioning uranium and 
is expelled through the nozzle to produce thrust. The moderator is 
used to reflect and reduce the energy of ( thermal ize) the neutrons 
which aids in sustaining the fission process. The pressure shell 
encloses the core, cavity, and moderator, while the radiator is 
used to remove any excess heat that builds up in the engine which 
cannot be removed by the hydrogen alone. 

Once a working knowledge of the engine had been acquired, it 
was then possible to compute various important engine character- 
istics. These computations were done using a computer program 
written during this study, which can be found listed in Appendix B. 
This program was used to generate the important parameters of the 
SRCGNR from general equations and equations given by Hr. Ragsdale 

•j 

for estimating the mass of an engine . The results obtained from 
this program may be found tabulated in Table 3.2. The program has 
four inputs that need to be determined before the other parameters 
may be calculated. These inputs are I $p , Fj , cavity diameter, D c , 
and the ratio of the volume of the gaseous uranium core to the 
volume of the cavity, V u /V c . For the data in Table 3.2, these 
values were selected to be I sp = 5000 seconds, F^ = 4.4 x 10^ 
newtons, D c = 3.3 meters, and V u /V c = 0.23. These figures were 
considered to be reasonable estimates of the operating conditions 
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Table 3.2 - Engine Parameters of a SRGCNR, = 4.4 x 10 3 N, 
SPECIFIED VALUES: 


I =5000.00 s Thrust=440000 . 00 
CALCULATED ENGINE PARAMETERS: 

ENGINE CAVITY CONDITIONS: 
T(H^) = 25165.49 K 
Cavity Pressure=981 . 25 atm 
Cavity Volume=18.82 nr 
Uranium Diameter=2 . 02 m 
Hydrogen Molecular Mass=0.759 


MASS FLOW RATES: 
H 2 -To-U 233 Ratio=393 . 89 
^=8.951 kg/s 


ENGINE DIMENSIONS: 

Pressure Shell Thickness=0 . 087 
Pressure Shell Mass=52597 kg 
Moderator Thickness=0 . 76 m 
Turbo Pump Mass=773 kg 
Radiator Mass=138463 kg 


VARIOUS ENGINE PARAMETERS: 
Reactor Power=16185 . 08 MW 
Specific Mass=0. 02224 kg/KW 
Thrust to Weight Ratio=0.188 


THROAT CONDITIONS: 
Temperature=22258 . 87 K 
Hydrogen Molecular Mass=0.843 
Area=0.0022 m 2 


EXIT CONDITIONS: 
Temperature=4120 . 51 K 
Hydrogen Molecular Mass=2.004 
Area=0 . 6489 m 2 

Exit-To-Throat Area Ratio=300 


N V a /V c =.23 D c =3.30 m 


T(U 233 ) = 76183 . 45 K 
Critical Mass U 233 = 40.12 Kg 
Uranium Volume=4.33 m 3 
Uranium Density=9.27 Kg/m 3 
Gamma =1 . 261 


Propel lant=8 . 97 kg/s 
U 233 =0.023 kg/s 


m 

Moderator Mass=46632 kg 
Exhaust Nozzle Mass=233 kg 
Total Engine Mass=238682 kg 


Radiated Power=971.10 MW 
Jet Power=10787 . 32 MW 
Engine Ef f iciency=0 . 6665 


Pressure=524 . 48 atm 
Gamma=l .246 
Velocity= 16542.23 m/ s 


Pressure=0 . 115 atm 
Gamma= 1 . 226 

Velocity= 49033.25 m/ s 
Mach Number= 6.50 
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of the SRGCNR engine 7 . With these inputs, the program was able to 
calculate the conditions in the cavity, in the throat, and at the 
exit of the engine. The formulas that were used to calculate these 
conditions are included with the program in Appendix B. Some of 
the more interesting results noted from Table 3.2 are a throat 
temperature of 22,260 K, a jet power of 10.8 gigawatts, a reactor 
power of 16.1 gigawatts, a thrust-to-weight ratio of 0.19, and an 
exit velocity of 49 km/s. This value of exit velocity is over ten 
times what a standard chemical rocket can attain today. The over- 
all mass of the engine was found to be 238.7 metric tonnes. 

Once the SRCGNR characteristics and values were obtained, it 
was then possible to proceed with other calculations that were 
related to the propulsion system. For example, the total number of 
engines, n, that would be required for any given mission had to be 
determined. This formula, given by 

n - T ?P F. (e 7 ^ - i) (3.3) 

^pr 

was found by replacing J in Equation 3.2 with nF^t pr . Once the 
number of engines was determined, it was possible to calculate the 
total engine mass from 

7 n„ - n m vl . (3.4) 

Equation 3.3 is plotted in Figure 3.4, which shows the relation- 
ship between the number of engines required and the mission AV for 
a range of different powered-flight times. This figure was 
generated using an I s? of 5000 seconds, an Fj of 4.4 x lO"* newtons. 
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I„ = 5000 s = 1 x 10* kg F = 4.4 x 10 5M 



Figure 3.4 - Number of Engines vs. AV with Powered-Flight Time 
a Parameter 


and an assumed dry mass of 1 x 10 8 kg. By using Equation 3.3, the 
number of engines for any selected mission parameters (Av, t pr , m,j rp ) 
may be determined. Once the total engine mass has been found, it 
was then possible to calculate the payload mass in the form of the 
payload mass ratio, given by 

Ay _ Ay 

JOl - e~ ^ - W 9 ( 1 - e" x -*). (3.5) 
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Figure 3.5 - Payload Mass Ratio vs. Engine Specific Impulse 
(I spe ) with AV as a Parameter. 
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tends to level out no matter how much the engine specific impulse 
is increased. Thus, an increase in the I jpe is only beneficial when 
operating in the region where the payload mass ratio substantially 
increases. There is little to be gained by increasing the thrust 
or the powered-flight time or decreasing the engine mass once the 
I Jpe has reached the point where the payload mass ratio begins to 
level off. Also, it may be noted from Equation 3.5 that there is 
a value of I $pe below which, the engine mass exceeds the entire 
payload mass, and thus the payload mass ratio will be negative. 

3.3 CONCLUSIONS 

The end result of the analysis presented in this chapter was 
to find equations that would yield the total mass of the propulsion 
system and propellant required. By using the equations and figures 
given here, these masses can be found for different mission specif- 
ications. This is a crucial step in the design of the Emerald City 
for specific missions. Two missions were designed for the Emerald 
City, using the equations derived here.. These missions can be 
found in Chapter 6. 

The preliminary propulsion system that was chosen was the 
space radiated open-cycle gas-core nuclear rocket engine. The 
benefits of this system are its high I sp and thrust values that can 
be obtained. However, this engine is still very conceptual in 
nature, and so it is not currently known if the idea will 
eventually lead to a functioning rocket engine or not. It can be 
seen from Table 3.2 that the temperature in the throat of the 
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nozzle is predicted to be around 22,000 K!. Clearly, advances will 
have to be made in the development of materials capable of 
withstanding this kind of thermal load. Another area in need of 
development is the confinement of the gaseous core. Currently, it 
is envisioned that the core will be confined by a vortex of inert 
gas. This area will need to be thoroughly researched if the 
concept of an open-cycle gas-core reactor is ever to reach 
fruition. 

There were two detrimental side effects noted by using the 
open-cycle engine, both of which are related. The first of these 
effects was the loss of uranium out of the cavity and through the 
nozzle. It was estimated that approximately 2 metric tonnes of 
uranium will be lost per engine per day of powered-flight time. 
The second side effect is the environmental and human factors 
concern of the radiation generated by this type of engine. Large 
amounts of radiative material are present in the exhaust plume 
which will result in massive shielding requirements for the crew. 
Both of these problems can be lessened by going to a light-bulb 
gas-core engine, which is a closed-cycle type of engine. However, 
this engine is not foreseen to have an I sp as high as that of the 
open cycle engine. The light bulb engine could prove useful if the 
scope of Project WISH were to ever be reduced, i.e. fewer people, 
less massive ship, smaller Av's... Until then, however, the open 
cycle gas-core engine will remain as the only type of nuclear 
engine capable of generating the high impulse-high I jp combination 
needed for Project WISH in the mid-21 s ^ century. 
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CHAPTER 4 


ATTITUDE CONTROL 


4.0 INTRODUCTION 

Attitude control is an important facet of Project WISH due to 
the fact that it will be necessary to damp out oscillations in 
attitude due to disturbances. Because the Emerald City will be 
spinning to produce artificial gravity (see Chapter 5), the system 
will be gyroscopic in nature. Therefore, stable equilibrium 
configurations had to be determined, which was done during the 
Phase I study conducted last year. Once these configurations were 
known, attitude control studies were conducted in order to find the 
control requirements needed to return the Emerald City to 
equilibrium following a disturbance. The methodologies used to 
perform these attitude control studies are the focus of this 
chapter . 

4.1 ATTITUDE DYNAMICS BACKGROUND 

For a body that is as large and complex as the Emerald City, 
there will be some structural deformations that arise due to the 
loads encountered. However, for the purpose of this preliminary 
investigation, the Emerald City was considered to be a rigid body, 
i.e. there are no deformations present. Although structural 
dynamics was not considered in this analysis, it was examined last 
year using a lumped-mass modelling method, and will again be 
considered in detail during the third year (1991-1992) of the 
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project. Since the goal of the attitude control study was to 
obtain an order of magnitude estimation of the control thruster 
requirements, it was felt that the rigid body assumption was 
adequate . 

As stated earlier, the stable equilibrium configurations for 

the Emerald City were determined from last year's vehicle dynamics 

2 

study . These stable configurations were found to be dependent on 
the parameters b and r in the following manner: 


1. 

b 

> 

1/r 

and 

b > 

4/r - 3 

2. 

b 

< 

1/r 

and 

r > 

1 

3. 

b 

< 

4/r 

- 3 

and 

r < 1. 


The parameters b and r are given by 

r = I z / I x = slenderness ratio 
b = n / Q = spin rate / orbital angular rate 
where I s and I x are the mass moments of inertia about the z- and x- 
axes (see schematic in Chapter 6), and n is the spin rate of the 
torus. Thus, it may be seen that the stability of the Emerald City 
depends upon three principal factors: mass moments of inertia 

(determined by ship configuration), spin rate (determined from 
artificial gravity considerations), and orbital angular velocity 
(determined by the nominal orbit). 

4.2 CONTROL SYSTEM REQUIREMENTS 

Once the stable configurations of the Emerald City are known, 
the control aspect of the problem may be examined. The main goal 
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of the attitude control study was to determine the values of 
thrust, root-mean-square power, and propellant mass required to 
damp out an initial disturbance. 

4.2.1 State Feedback Control Design 

The first step in achieving this goal was to determine the 
state-space equations that model the motion of the Emerald City. 
This equation was expressed in the non-dimensional state-space 
formulation 8 

x - [A] [ x] + [B][f] f 4 * 1 ) 

where 


[A] 


0 0 

0 0 

(l-r) 0 

0 (l-r) 


1 0 

0 1 

0 (2-r) 

(r-2) 0 


[B] 


0 O' 
0 0 

1 0 

p 1. 


/S * * *■ 

W3 - [0 t 0 2 0! 0 2 J 


(4.2) 


and [T] is the non-dimensional torque provided by the control 
thrusters ( ~ denotes non-dimensional ized quantities). Equation 4.1 
is used to simulate the state response of the Emerald City for 
given ship configurations, thrust input configurations, and initial 
disturbances. A linear quadratic regulator control program 
(MATLAB 8 software) was used to obtain the non-dimensional state 




response that minimizes the control design performance index, 
defined by 

CDPI - -if" (£ 7 J? + f T f )dv. (4.3) 

2 Jo 

By minimizing this performance index, the program minimizes the 
total control effort and the kinetic energy of the system and 
yields a linear state feedback control law 

[T] - [G][x] (4-4) 

where G is the control feedback gain matrix minimizing the CDPI in 
Equation 4.3. 

4.2.2 Attitude Control Power Required 

For the control torques of Equation 4.4, the non-dimensional 
control power consumed can be found by the non-dimensional power 
S*, 

5- - rf T r T f- 1 fdx (4.5) 

Jo 

which yields 

5* - XjP m X 0 (4.6) 

for any initial disturbance state x fl . The power matrix P* is the 

solution of the linear matrix equation 

a/ 2 P * + P'Acj + G fr T - 0 (4.7) 
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where 


and the non-dimensional 


Ac/ “ 
thruster 


A + BG (4.8) 

distribution matrix is defined by 


$ - Dl 1 [D] (4.9) 

in which is the torus diameter and the elements of [D] are the 
moment arms of respective input forces yielding moments about the 
centroidal body axes x and y of the vehicle which are parallel to 
the torus plane. The required quantities and matricies in 

Equations 4.4 through 4.9 can be evaluated by using any suitable 
control software such as the MATLAB Control Toolbox mentioned 
previously . 

The various dimensional and non-dimensional quantities are 
related by 


[ »i e 2 8, 8 2 ] - [ §, S 2 e, e 2 ] 


10 0 0 
0 10 0 
0 0 n 0 
o o o n 


(4.10) 


(4.11) 


p T r n2 

F 1 F i n 


(4.12) 


T - f I, n 2 . 


(4.13) 
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The non-dimensional control thrust forces can be obtained via 


p - fr'f ( 4 . 14 ) 

In all of the above equations, a ~ above a variable denotes non- 
dimensional quantities. The root-mean-square power required to 
damp out a disturbance can then be found from the value of S* , 
computed via Equations 4.6 and 4.7. The root-mean-square control 
power required is then given by the formula 




( 4 . 15 ) 


where v gI is the exit velocity of the control thrusters. 


4.2.3 Attitude Control Propellant Requirement 

The propellant mass required for the control thrusters can be 
found by the formula 

where f is the total number of inputs used. The integral in 
Equation 4.16 is found by numerically integrating the non- 
dimensional thrust input versus non-dimensional time plots (see for 
example Figure 6.4 for specific mission designs) 


4.2.4 Attitude Control Thruster Configuration 

For a body with the size of the Emerald City, one would need 
hundreds, perhaps thousands, of control thrusters to exert attitude 
control torques about the body axes. The number would depend on 
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the thrust output capabilities of the control thrusters. Taking 
into account so many thrusters individually in the attitude control 
problem study would result in a very large D (thruster distrib- 
ution) matrix. However, because the Emerald City is assumed rigid, 
one can instead consider only a few clusters of thrusters, where 
each cluster is represented as if it is a single thrust force 
obtained as the resultant force of hundreds of thrusters each 
producing identical thrust-time histories in a given cluster 
location. Thus, each cluster location acts as one input in the 
control study. After finding the total thrust-time history at a 
particular location through the control study, one can then 
apportion equally the required thrust to individual thrusters. 
Thus, for given thrust capabilities of thrusters, the total number 
of thrusters required in a particular cluster can be found 
subsequently. For example, in Figure 6.5, only three cluster 
locations are shown as if there are only three control thrust 
inputs, F^ , Fj, and F 3 , and the attitude control requirements are 
computed for that particular input configuration. How many 
individual thrusters may exist at each cluster becomes immaterial 
at this point and the control study can proceed. 


4.3 ATTITUDE CONTROL DESIGN METHODOLOGY 

The methodology developed in this chapter can be used to study 
the attitude control of the Emerald City for a wide range of ship 
configurations, thruster configurations, and initial disturbances. 
Two particular scenarios were examined and the control system 
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requirements were determined using the methods provided above. The 
results for particular mission studies can be found in Chapter 6. 
The outline of the methodology as described above is reviewed here 
in step-by-step form. 

1. Choose a cluster configuration, characterized by Equation 
4.9, to be studied and determine the slenderness ratio, r, 
from the ship configuration. Select an initial mean 
disturbance to investigate. 

2. Use MATLAB (or a similar linear quadratic control 
regulator program) to determine the feedback gain, G, and 
the state response, x, to Equation 4.1 that minimizes the 
control design performance measure of Equation 4.3. 

3. Use MATLAB to find the control torque, T, and the 
resultant cluster inputs, F (via Equation 4.14), versus 
time profiles, the non-dimensional power parameter, S* , 
and observe the non-dimensional control time, t c , from the 
x response plots. 

4. The resultant cluster thrusts required, F' , for the 
control system can be found from F and dimensional izing 
according to Equation 4.12. The root-mean square power 
and propellant required can be found from Equations 4.15 
and 4.16, respectively. 

One of the more important results obtained from the attitude 
control study was that the root-mean-square power required to damp 
out a disturbance is proportional to the square of the spin rate. 
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as shown in Equation 4.15. However, the torus diameter, D^, can be 
written in terms of the spin rate and the g-level to be provided 
(see Chapter 5) as 


Dt 


n 9 g [ 

60 V 

2nl ( 

2%) 


(4.17) 


where the spin rate n s is now in revolutions per minute (rpm's). 
Inserting this equation into Equation 4.7 and converting the spin 
rate to rpm's yields 


1 2*1 

a l z n s 

S* 

l 60 j 

1 a s9 > 

T c 


(4.18) 


Equation 4.18 clearly shows the benefit of keeping the spin rate 
below 1 rpm (since for n < 1, n* < n) in order to decrease the 
power requirements for attitude control. This result was taken 
into account in human factors studies when the torus was designed 
(see Chapter 5 ) . 


4.4 CONCLUSION 

This year, the work on the dynamics and control aspect of 
Project WISH focused on determining the attitude control system 
requirements. A methodology was established to perform attitude 
control studies for a wide variety of missions. This methodology, 
which is summarized above, was used to design the attitude control 
system for the two missions of the Emerald City. The results of 
these mission designs can be found in Chapter 6. 


4.9 


CHAPTER 5 


HUMAN FACTORS 


5.0 INTRODUCTION 

Maintaining a habitable environment in space presents problems 
with varying degrees of difficulty depending upon its size and how 
long that environment is to be maintained. Something such as food 
supply, which is a relatively minor concern on a shuttle mission, 
becomes a problem that could dominate the design of a spacecraft 
with a crew of between 500 and 1000 and an operational lifetime of 
fifty years. Because of the magnitude of the Emerald City, human 
factors will play an important role in its design. This chapter 
covers three topics that relate directly to the design of the 
spacecraft: the life support system, the need for artificial 
gravity, and the problem of radiation. 

5.1 LIFE SUPPORT SYSTEM 

Last year, the different functions of life support were 
approached using mechanical systems^. Separate systems were used 
for air revitalization, temperature and humidity control, waste 
management, and so forth. However, there are two major faults with 
using mechanical systems. Because of the three year plus mission 
times expected for Emerald City, all of these systems would have to 
run (unrealistically) with one-hundred percent efficiency. This is 
necessary to prevent compounds such as air or water from being 
recombined into products that cannot be recycled. More importantly 
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though, is that since these systems are mechanical, they are 
subject to breakdown and will eventual ly need spare parts. Storing 
these parts or making provisions for their manufacture would waste 
valuable space and add extra mass to the ship. For these reasons, 
another concept was explored for the primary life support system: 
biospheres . 

A biosphere is a totally enclosed ecological system, with 
energy as its only input. The collective biomass of the earth is 
referred to as Biosphere I. Since all air, water, and waste is 
recycled internally, life can be sustained indefinitely. Clair 
Folsome created the first human-made biosphere inside a lab 
flask*®. The flask contains only simple organisms, but it has been 
alive since it was sealed in 1968. What is proposed for Project 
WISH is that a scaled and more efficient version of Biosphere I be 
used as the primary life support system. Research is already being 
done to create biospheres large enough for human occupation. The 
Biosphere II project*® in Oracle, Arizona has constructed an 
environmental facility sealed off from the atmosphere that uses its 
land and water four and ten times as efficiently as Biosphere I**. 
Total recycling of all atmospheric carbon dioxide takes place in 
half a day instead of eight years. The 150,000 cubic meter 
facility is designed as an environment for eight people who were 
recently sealed inside for a two year study. This figures to 
roughly 19,000 cubic meters per person. Even though this number 
would probably shrink if research were devoted to reducing the 
volume needed per person, it was used for sizing the crew section 
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to allow extra space for manufacturing, food processing, or other 
needs that have not yet been considered. The equation for the 
volume required in the torus section was then given by 

^eoru* ” 19000W (5.1) 

where N is the number of crew members on the Emerald City. 
Research is currently aimed at understanding more about Biosphere 
I and to possibly extend man made biospheres onto other planets. 
It was believed that the concept of biospheres could legitimately 
be considered for a spacecraft with the life expectancy and 
magnitude of Project WISH. 

5.2 ARTIFICIAL GRAVITY AND SHIP DESIGN 

Of the little information known about long term effects of 
zero-g on the body, none of it is very positive. Without gravity, 
bones decalcify and muscles weaken and atrophy. A condition that 
could be described as "space anemia" also occurs . This is due 
to a drop in the production of red blood cells. This otherwise 
harmless condition could pose a serious health problem if an 
astronaut were injured or became ill. To counter these effects, 
the crew section of Emerald City was designed as a torus rotating 
about a central hub to provide some level of artificial gravity. 
Even though providing 1 g would eliminate these effects, consider- 
able power and mass savings can be made if an acceptable lower g- 
limit can be found. The relationship between rotation rate, 
distance from the spinning axis, and g-level is shown in Figure 
5.1. Since the volume requirement (from Section 5.1) and shape for 
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the crew section is known, the actual geometry is determined from 
the g-level and spin rate desired. Several possible geometries are 
listed in Table 5.1. Selection of the major radius of the torus, 
R (see Figure 6.1), is dependant on the constraints set for g-level 
and spin rate. One revolution per minute is used as an upper limit 
on spin rate because the power needed to maintain the attitude of 
the Emerald City is directly proportional to spin rate to the 
fourth power. This is discussed in greater detail in Chapter 4. 
Another reason behind this limit is human endurance. Most humans 
cannot endure extended periods of rotation at rates higher than 1 
rpm 1 ^. Although there is no universal agreement on this limit, 1 
rpm was chosen as it was the most conservative limit found. To 
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Table 5.1 - Possible Torus Geometries 
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540 
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11.1 
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1.12 
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580 

29.0 

10.0 

609 

1.08 

1.01 

0.94 

620 

28.0 

9.0 

648 

1.05 

0.98 

0.91 

reduce the mass 

660 

27.2 

8.2 

687 

1.02 

0.95 

0.88 

700 

26.4 

7.5 

726 

0.99 

0.93 

0.86 


further, but more 


information on the effects of long term exposure to lower gravity 
levels is needed before a lower g-limit can be set. 

Another factor that contributes to the mass of the torus is 
the fact that it must hold an atmosphere. If a stressed skin is 
used, the necessary skin thickness and structural mass (M ss ) of a 
torus, with major radius R and minor radius r, is found using the 
formulas shown on the top of the next page^ . The first two 
formulas find the skin thicknesses needed to withstand the prin- 
cipal stresses on the surface of the torus. Generally, t^ oop is 
greater than tneridional ' so fc hoop ^ s use d to find the structural mass. 
Structural mass is calculated for several configurations later in 
this chapter using the working stress ( a w ) and density (p) of 
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2 R % D 

°„-p* R 
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aluminum. P g is the pressure due to the load distribution inside 
the torus. One way to reduce the structural mass is to lower the 
atmospheric pressure inside. Doing this would also result in 
effects similar to living at higher altitudes on Earth: the 
coughing mechanism would be less effective, a normal speaking voice 
would not carry as far, boiling temperature would be lowered, etc. 
Since structural mass is relatively minor compared to the whole 
ship, a pressure of one atmosphere is used in the crew section for 
all calculations. 


5.3 RADIATION AND SHIELDING 

The principle risk of living in space is from radiation 
induced cancer. The National Council on Radiation Protection and 
Measurements^ has recommended a career limit for whole body 
exposure with a lifetime three percent excess risk of fatal cancer. 
Figure 5.2 is based on age at the start of exposure with a ten year 
career assumed. This recommendation, however, is made not for 
exploratory missions, such as a Mars mission or those that would be 
undertaken by the Emerald City, but for multiple missions of 
approximately 90 days over a 10 year period. What is needed for 


Project WISH is a limit 
for constant exposure. 
Table 5.2 shows the 
limits recommended for 
constant public and 
occupational exposure in 
the United States. The 
occupational annual 
exposure limit of 5 rem 
per year would be the 
limit used on the Emer- 
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Figure 5.2 - Career Limit vs. Age 


major sources of radiation for Project WISH: galactic cosmic 
radiation, radiation from solar flare events, the propulsion 
system, and the power system. Of these four, only the last two can 
be reduced by moving the source away from the crew or by covering 
the source with a shield. Cosmic rays are present throughout the 


Table 5.2 - Continuous Dose Limits 


Occupational 

Annua 1 5 . 0 rem 

Lifetime (guideline) Age*1.0 rem 

Planned or special 10.0 rem 

emergency 

Pub lie 

Annual , continuous 0 . 1 rem 

Annual, occasional 0.5 rem 


solar system at 
relatively constant 
intensities, decreasing 
only during solar flare 


activity 


Unshielded 


dose rates from cosmic 
rays are estimated to be 


from 30 to 50 rem per 
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year*®, while large solar flares are capable of delivering doses 
in excess of 1000 rem in only a week*®. Since shielding these 
sources is impossible, the only way to protect from them is to 
shield the crew section. Radiation from the power system and the 
gaseous uranium core of the propulsion system can be attenuated by 
separating them from the crew and by partial shielding (see Figure 
5.3). Unfortunately, the expanding plume of exhaust from the gas- 



core rocket contains fission fragments that emit gamma radiation 
back towards the crew section. Increasing the distance between the 
crew and the nozzle exit will reduce the dose rate, but structural 
limits will only allow so much separation. The only other means of 
attenuation is to shield the crew section. The actual dose rate 
was calculated with the parameters of the gas-core engine from 
Chapter 3 using the formulas given by Charles Masser . It is 
primarily dependant on the amount of separation between the crew 
and the engine, and the amount of time that fission fragments 
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Pigure 5.4 - Dose Rate vs. Separation for Varying Fission 
Fragment Retention Times 

remain in the core before exiting the nozzle. Figure 5.4 shows 
this relationship for three different fission fragment retention 
times. The only estimate for the fission fragment retention time 
of the gas-core engine comes from Masser, who estimates it to be 
about 100 seconds. Note that the numbers in Figure 5.4 are only 
for one engine. It has been assumed that for the case of multiple 
engines, the dose rate will be additive. Although this assumption 
is probably not accurate, it is one that had to be made because of 
the limited knowledge of radiation by the members of this project. 

Galactic cosmic radiation is considered to be the primary 
source of radiation because it is always present and cannot be 
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shielded at the source. No matter what their intensity, solar 
flares and engine exhaust will only be present for periods of about 
a week or two at a time. Because of this, the selection of a 
primary form of shielding was based on protecting against galactic 
cosmic radiation. The most common method of shielding is passive 
shielding. Passive shielding, known to be effective, simply puts 
a physical barrier between the source and the target. A more 
theoretical type of shielding is active shielding. It uses either 
an electromagnetic or electrostatic field to repel incoming 
radiation. A drawback common to both field types is that, 
depending on how much power they receive, they only shield 
radiations up to a certain energy level. Electrostatic fields were 
eliminated because a potential on the order of 3 GV is needed on 
the surface of the torus . This potential would limit extra- 
vehicular activity and also generate lethal bremsstrahlung 
radiation fields inside the torus. An electromagnetic shield was 
rejected because higher energy radiation could become trapped 
inside the field. This would create something similar to the Van 
Allen belts around the Earth only much more dangerous considering 
the smaller distances involved. Another common failing of the 
shields is that they only stop charged particles. Gamma rays and 
neutrons, the more destructive types of radiation, are not 
affected . 

With only passive shielding remaining, materials that could 
best attenuate radiation and still save mass were explored. Town- 
shend compares different materials and finds liquid hydrogen to be 
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very effective**. Liquid hydrogen is more effective than other 
materials because it contains no neutrons for the incident 
radiation to scatter upon collision. Table 5.3 shows that 100 
grams per square centimeter will reduce the dose to below 5 rem per 
year. One hundred grams per square centimeter of liquid hydrogen 
corresponds to a shield thickness of 14 meters. The mass of this 
type of shield, as well as the structural mass, is found for four 
different torus configurations in Table 5.4. Notice that even if 
the crew size is cut in half, the masses for the torus do not 
significantly drop. This is because of the geometrical nature of 
a torus, a large volume change results in a relatively minor change 
in surface area. 

5.4 CONCLUSION 

Even though all the considerations in this chapter have not 
been solved, enough has been determined to create a pattern that 
can be used to make an initial design for the Emerald City. The 
results of research in the areas of biospherics, microgravity, 
radiation, and other fields not yet considered will set the limits 
that help finalize any design. The information presented in this 
chapter will be used for a preliminary design of the Emerald City 
for two particular missions, which may be found in Chapter 6. 
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Table 5.3 - Solar Minimum Galactic Cosmic Ray Depth Dose 
Equivalent in Tissue as a function of Particle Type and LHj Shield 
Thickness 
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Table 5.4 - Structural and Shield Mass Estimates 


crew 

size 

major 

radius 

minor 

radius 

surface 

area 

shield 

mass 

structural 

mass 

skin 

thickness 


(m) 

(m) 

" U) 

(k«) 

(k«) 

(cm) 

1000 

500 

44.2 

8 . 725e+05 

8.725e+08 

4. 101e+07 

1.741 

1000 

700 

37.3 

1.031e+06 

1.031e+09 

5.052e+07 

1.815 

500 

520 

30.6 

6.282e+05 

6.282e+08 

2.407e+07 

1.419 

500 

700 

26.4 

7.296e+05 

7 . 296e+08 

3.027e+07 

1.537 
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CHAPTER 6 


REPRESENTATIVE MISSION DESIGN 


6.0 INTRODUCTION 

Once the methodologies of the previous chapters have been 
formulated, two sample missions were examined in order to obtain 
some preliminary sizing and mass estimates for the Emerald City. 
A Saturn Envelope mission and the currently popular Earth-to-Mars 
mission were chosen as representative scenarios for this 
preliminary design. The design for these missions were determined 
by incorporating the orbital mechanics, propulsion, attitude 
control and human factors subsystems. Orbital mechanics (Chapter 
2) was used to determine the Av's required for the selected 
missions. The propulsion system study (Chapter 3) was used to 
determine the number of engines and the total engine and propellant 
masses. These were then used to size the Emerald City in terms of 
volume of propellant, engine spacing and location, length of 
central pole, etc. Human factors (Chapter 5) was used to determine 
the sizing of the torus, and the attitude control study (Chapter 4) 
was used to calculate the control thruster requirements. With all 
of the work laid down in the previous chapters in mind, it is now 
possible to proceed with the design of both missions. 

6.1 PROCEDURE 

The procedure that was followed in the design of these 
missions is listed sequentially as follows. 
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1. Obtain the AV required for each mission from orbital 
mechanics . 

2. Set values for the specific impulse, I jp , and the thrust 
per engine, F^, for the propulsion system. The mass of 
each engine, m^, is then found from the computer program 
of Appendix B. 

3. Determine the torus mass, n»^ oras , from the masses of the 
biosphere, the shield, and the torus structure. These 
masses are dependant upon the number of crew members 
selected. The dry mass, m^, is then estimated by making 
allowances for the engine, cargo, and remaining subsystem 
masses . 

4. Calculate the number of engines required and the total 
engine mass, m 4 , for each mission. 

5. Calculate the payload mass, mj, for the missions. This 
mass is defined as m^ - m torflS - and includes cargo 

as well as other systems that have yet to be sized (heat 
rejection, power, tankage, etc.). 

6. Calculate the total initial mass, m 0 , and the propellant 

mass, mp. The propellant mass ratio, m^ / n^, and the 
payload mass ratio, mj / , are then found. 

7. Calculate the volume of propellant required for the 
mission, V p . Then, by estimating a stem radius, r p , 
the height of propellant, and hence the minimum stem 
length, , can be found. 
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Neptune and Pluto cannot be reached at all in flight times of up to 
five years with this AV ceiling. It should be noted that many of 
the transfers to Saturn can be made with AV's considerably less 
than 50 km/s round trip, but for this mission, it was assumed that 
the Emerald City was in a poor orientation (AV-wise) with Saturn 
and that a three year flight time was essential. A possible 
motivation for such a mission might be that the Emerald City is 
required to support the first manned exploration of Saturn's moon 
Titan, and that the planetary alignments and limitations on the 
interplanetary transfer ship require a particular launch date from 
the Martian colony and arrival time at Saturn. 

For this mission, the individual AV’s and the powered-flight 
times required for the four separate burns were not considered, but 
only the total Av and powered-flight time. The number of engines 
was found based on these total values and a given total powered- 
flight time of 20 days, but it should be kept in mind that the 
powered-flight times would not be the same for all four burns, even 
if the AV required would be the same. This is because the mass of 
the Emerald City is smaller for each of the burns, and hence the 
longest powered-flight time would be the first one. 

The results for the Saturn mission can be found in Tables 6.1. 
Some important values to note are the total initial mass, 4.16 x 

ft 

10 kilograms, the number of engines, 172, the torus radius, 700 
meters, and the height of the pole section, 1270 meters. The 
payload mass ratio for this mission was found to be 0.083, or 8.3% 
of the total initial mass is payload. It must be remembered, 
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however, that the payload mass as defined here is the mass of the 
cargo plus the mass for any other systems that were not considered 
during Phase II (power systems, thermal rejection systems, tankage 
mass , etc . ) , 

6.2.2 Earth-to-Mars Mission 

The second mission that was examined was a transfer from 1 AU 
to Mars. This mission represents a more common mission that the 
Emerald City might be required to undertake, perhaps carrying new 
colonists and supplies that had been shuttled to the Emerald City, 
just outside the Earth's sphere of influence, to the ever expanding 
Martian colonies. For this mission, the AV required was broken 
down into the departure and arrival AV's. The mission was 
scheduled to leave Earth on November 7, 2050 and the flight time 
was determined by minimizing the AV for this launch date using the 
MULIMP program. The AV ' s were found to be 5.1 km/ s and 7.5 km/ s 
for a total AV of 12.6 km/s, and the flight time was found to be 
1.52 years. The first powered-flight time was set at 5 days and 
the number of engines required was determined. Then, keeping this 
number of engines constant, the second powered- f 1 ight time was 
found. From Table 6.1, it can be seen that the results of this 
were found to be 33 engines, a 5 day initial powered-flight time, 
and a 6.53 day terminal powered- f 1 ight time. Some of the other 
important design values for this mission are the initial mass, 
1.295 x 10 5 kg, the torus radius, 700 meters, and the pole height, 
569 meters. 
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Table 6.1 


Summary of Design Variables for Sample Missions 



Saturn Mission 

Earth to Mars 

People 

1000 

500 

AV^ al (km/s) 

50 

12 . 6 

Av, (km/s) 

- 

5.1 

Av, (km/s) 

- 

7.5 

F ; (Newtons) 

4.4 x 10 5 

4.4 x 10 5 

I eB (seconds) 

5000 

5000 

(days) 

20 

11.53 

t nr , (days) 

- 

5 

t nr , (days) 

- 

6.53 

number of engines 

172 

33 

(kg) 

1.5 x 10 9 

1 x 10 9 

%«r„c (kg) 

1.111 x 10 9 

7.746 x 10 8 

m„ (kg) 

4.4376 x 10 7 

8.514 x 10 6 

m, ( kg ) 

3.457 x 10 8 

2.165 x 10 8 

m^ (kg) 

2.658 x 10 9 

2.947 x 10 8 

m„ (kg) 

4.16 x 10^ 

1.295 x 10^ 

m, / m„ 

0 . 083 

0.167 

m. / m n 

0 . 639 

0.228 

V !R , (m 8 ) 

3.742 x 10 7 

4.151 x 10 6 

(m) 

100 

50 

h n (m) 

1200 

528 

h (m) 

1270 

569 

R (m) 

700 

700 

(m) 

37 

26 

I, (kg m 2 ) 

5.6144 x 10 14 

3.8061 x 10 14 

I T , I ? (kg m 2 ) 

9.9568 x 10 14 

2.1731 x 10 14 

r (I, / I T ) 

0 . 563876 

1.75154 j 

max/min g-levels 

0.8/0.72 

0 . 8/0 .74 

n (spin rate, rpm) 

0 . 99 

0.99 
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Figure 6.1 - Schematic of the Emerald City 

6.2.3 Attitude Control System Requirements 

An attitude control system was designed for the two particular 
missions examined for a sample initial disturbance and a sample 
thruster configuration. Using the methodology developed in Chapter 
4, the root-mean-square power, propellant mass, maximum thrust, and 
control time was determined, as may be seen in Table 6.2. Also, 
the state response, control torque profile, and thrust profile for 


6.7 





the Saturn Envelope mission configuration are shown in Figures 6.2- 
6.4 for the sample thruster configuration shown in Figure 6.5. 


Table 6.2 - Attitude Control System Requirements 


I* / I* 

0.56 

1 .75 
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6.176 x lO 12 

7.86 x 10 U 

"V (kg) 

2.105 x lO 7 

3.458 x 10 6 

1 mav _ 

2.988 x lO 9 

1.911 x lO 9 

^ control ^ 3 \ 
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Figure 6.2 - Non-Dimensional State Response (r = 0.5 ) 
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Figure 6.3 - Control Torque Profile (r = 0.5 ) 
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Figure 6.5 - Sample Thruster Configuration 


6.3 CONCLUSION 

In this chapter, two particular missions were examined for the 
Emerald City. However, using the methodology established here, a 
wide variety of missions may be looked at and the configuration for 
the Emerald City may be determined. In the future, different 
mission parameters, such as crew size, Av requirements, sample 
initial disturbances, attitude thruster configurations, etc..., may 
be examined. It is felt that this methodology is the main 
accomplishment of the Phase II design of Project WISH. 
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CHAPTER 7 


CONCLUSION 


7.0 PROJECT SUMMARY 

Phase II of Project WISH was able to accomplish a great deal 
towards the design of the Emerald City. In orbital mechanics, a 
nominal orbit of 4 AU ' s was selected, and the operational envelope 
for the Emerald City was determined. An open-cycle gas-core 
nuclear thermal propulsion system was chosen as the main propulsion 
unit, and a detailed study was conducted on this type of engine. 
For attitude control, a methodology was established to determine 
the thruster requirements needed to damp out a disturbance 
experienced by the Emerald City. Human factors looked at the 
problems of life support, artificial gravity, and radiation 
shielding. Biospheres were examined as the primary method to 
sustain life, while the torus size requirements were determined for 
the different levels of gravity that the crew would need. 
Radiation shielding proved to be one of the major design problems 
encountered, and mass estimates were obtained for shielding from 
cosmic radiation. The ultimate result of the Phase II study was 
the generation of Tables 6.1 and 6.2. These tables incorporated 
information from all of the subsystems studied, and provide data on 
the estimated size and mass requirements of the Emerald City. By 
examining these figures, one can realize the magnitude of a ship 
such as the Emerald City. 



concessions to the unknown, we had best frankly recognize that 
irrational faith in the value of a new capability often plays an 
essential part in the early evolution of a new concept. Certainly, 
at that point, the lack of a recognized application is a demon- 
strably unsound basis for the concept's rejection." 

Table 7.1 - Envisioned Time Line for Project WISH 

1989 - President Bush announces Mars Initiative to reach the 
Red Planet in 30 years. 

1991 - NASA presents a program to send sensing probes 

throughout the solar system. Projects include the 
return to the inner solar system, exploration of the 
asteroid belt, and further missions to the outer 
planets . 

1998 - Space Station Freedom becomes operational . 

2001 - Heavy Lift Launch Vehicle makes its maiden flight. 

2004 - Construction begins on a near-geosynchronous Earth 
orbiting space station. 

2009 - U.S. returns to the moon. 

2011 - Construction begins on the Moon Base. 

- Maiden flight of National Aerospace Plane. 

2015 - Moon Base becomes fully operational. 

2019 - First manned mission to Mars. 

2023 - First living modules constructed on Mars. 

- Construction begins on Reusable Interplanetary Ships 
(R. I. S.'s) for carrying personnel and cargo. 

2028 - Mars Base becomes fully operational. 

2040 - Implementation of Project WISH. 

2050 - The Emerald City becomes operational. 
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APPENDIX A 


EQUATIONS FOR AV MINIMIZATION PROBLEM 


The objective function given for the Av minimization problem 
set up in Chapter 2 was given by Equation 2.4 

AV t - A V t + AV 2 - V x - V a + V e - V 2 (2.4) 

where AVj and AVj are the velocity changes along the thrust vector 
direction, which will be kept constant during each powered-flight 
phase. The constraints were given in Chapter 2 as 

1 * t f£ + Sri + Sr2 < T0F »ax 

2 - a o < a o,«ai 

3 - tpri + < t pr,aai 

4 ’ x f,EC = x f , TO * 

a o' fc prl' fc pr2' > °‘ 

Thus, the problem is to determine the relationships between these 
values and the design variables of the problem. 

Since a circular nominal orbit and target planet orbit has 
been assumed, Vg and v^ are given by the formulas for circular 
orbital velocity, 

( A. 1 ) 




(A. 2) 
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which are known for a given nominal orbit and target planet. 
However, finding the velocities , the velocity at the end of the 
first powered- f 1 ight phase, and Vj, the velocity at the beginning 
of the second powered-flight phase, is more involved. To find the 
equations for these velocities, the kinematic relationship 

dv - a-dt ( A. 3 ) 

was integrated along the direction of thrust, with an acceleration 
term in this direction given by 


a( t) 



(A. 4) 


This term arises from the fact that as the Emerald City is burning 
with a constant thrust, it is expending mass, and thus the 
acceleration will increase with time according to Equation A. 4. 
Inserting Equation A. 4 into Equation A. 3 and performing the 
integration yields 


v 


v i ~ I at>sr 1^(1 - 


^5e). 




(A. 5) 


which becomes 

Ai; - v„ - - - 1^, InU - <A- 

for the first powered-flight portion of the transfer. Similarly, 
for the second powered-flight portion of the transfer, the velocity 
change is given by 
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A V V„, - v„ - -l^g la<l - -2a^) 

where (».7) 


The subscript F in Equations A. 6 and A. 7 denotes the components of 
the respective velocity vectors projected along the thrust 
direction (it is assumed that the thrust angles, P 12 , wil1 not be 
changed during a powered-flight phase). It can easily be seen that 
by using Equations A. 6 and A. 7, the objective function of Equation 
2.4 can be written as 

Av r - -r^tlnd - -SgBtB)+ln(i - H . (a.8) 


where a g2 is not an independent variable, but given in Equation A. 7. 

The application of the constraint functions needed many more 
equations and relationships. First, continuing in the same manner 
as was used to find Equation A. 5, the relationship 

dx ” v'dt ( A. 9 ) 

was integrated with the results from Equation A. 5 inserted for v, 
and the result was found to be 


c - x -*■ v t 
C l,f X o. 2 + V 0,2 c prl,2 + 


* 0,02 


where 


(A. 10) 


4 


- (1 - ) [ln(l - ( a °.° 2 C pn ,2 ) _ l] + 1 . 

1,2 -Tsp 3 Imp 3 
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Equation A. 10 yields the position along the direction of the thrust 
vector and the subscript l,f means at position 1 or f, etc. Thus, 
armed with Equations A. 5, A. 10, and some fundamental relationships 
from astrodynamics , the constraints can be related to the design 
variables of the problem in the following manner. First, the 
inertial X- and Y-components of the velocity and radius vectors at 
the end of the first powered-flight phase can be found by applying 
coordinate rotation matrices to Equations A. 5 and A. 10. These 
equations are given by 




a„fc. 


-Jr sin0 o - - e o )ln(l 


v iy “ 


-Jr cose o - r ap grcos(P 1 - 6 0 )ln(l - ) 


1 


(A. 11) 


\ " five + 


and 


r ix “ r 0 C0Be 0 - 


-£■ t^sinS „ * B in(8 1 - 8„>« 


N 


r iy “ Xo 3in0 o + 


t_ lC ose o + • (J qP sr) . a coa<p a - 0 O )* 

d _ 


jxrl 


(A. 12) 


* - (1 - a ° t P rl ) [ln(l - a °^ vzl ) - 1] + 1 

I sp& 


-1 - 


lx 


+ I- 


i y 


where the known variables are r Q and I sp and the design variables 
are 0 O , t prl , a Q , and (3^. Once the velocity and position vectors at 
the initial point of the free-flight trajectory (point 1) are 
known, the transfer orbit is defined by the following equations: 
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(A. 13) 


e *~ V [(V 1 " ^ r i* ~ (r ix Vix + r xy V ly > V lxl 

e y - -J- [ ( vi 2 - r iy - (r^ + r ly v ly ) v ly ] 

e - \jel + e y 

Y - cos' 1 ( — ) ( quadrant ) (A. 14) 

e 

Z - iL - JL (A. 15) 

2 r x 

3 “ _ 2l (A. 16) 

Y t - cos'M— — — ) (quadrant) (A. 17) 

3 JO* 


h - <Ja (1 - e 2 ) |i 


(A. 18) 


In the above equations, e is the eccentricity of the transfer 
orbit, x Q is the direction of the eccentricity vector, E is the 
specific energy of the transfer orbit, a is the semi -major axis of 
the transfer orbit, tj is the true anomaly at the end of the first 
powered-flight time, and h is the angular momentum of the transfer 
orbit. Once the transfer orbit is defined through these equations, 
the variables at the terminal free-flight point (point 2) may be 
calculated. The design variable that was designated was y^ , which 
is the true anomaly at point 2. This then defines the free-flight 
time of the transfer according to 
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(A. 19) 


tft ~ 


— (2*le + (E~ - e sinE,) - (E^ - e sini^) ] 


where E is the eccentric anomaly, defined by 


. , e + co8Y n , 

E n - cos' 1 (- ) 

n 1 + e cosy,, 


(A. 20) 


and k is the number of times the periapsis is passed during the 
transfer. The rest of the values at point 2 are given by 




a(l - e 2 ) 

1 + e cosy 3 


(A. 21) 


*2 - 


N 


2(E + -^) 


(A. 22) 


0 - [ 1 1 1 o Ya - INTEGER ( 7- ) ] •360 < 

2 360 


(A. 23) 


360" 


<b, - cos' 1 ( — — ) (quadrant) 

T * r rr 


-^2 


(A. 24) 


- r a cos0 2 
r 2y - r 2 sind 2 


(A. 25) 


v 2x ■ -v 2 sin(0 3 - + 2 ) 
v 2y “ < e 2 - * 2 > • 


(A. 26) 
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The equations governing the final powered-flight phase of the 
trajectory are found by applying coordinate rotation matrices to 
Equations A. 5 and A. 10, similar to the way that Equations A. 11 and 
A. 12 were obtained. These equations were found to be 


a„,fc r 


v fx~ -v 2 sin(0 2 - $ 2 ) - r^gsin ((J 2 - 0 2 )ln(l * 

Vfy ~ v 2 cos (0 2 - <j> 2 ) - Igpgcos (P 2 - 0 2 )ln(l ‘ijg* * 


(A. 27) 


and 

r fx - r 2 co80 2 - v 3 t pr2 ain(B 2 - <fr 2 ) + “ ® 2)®2 

Tfy - r 2 sin0 2 + V 2 t^cos (0 2 - ^ 2 ) + - 008 ( Pa " 0 2 > *2 ( A . 2 8 ) 

a - (1 - fg&Si ) [ln(l - - 1] + 1. 

2 r «pSr 

The above equations will play an important part in the 

determination of the constraints on the total time-of -f 1 iqht and 

the rendezvous requirement. For the time-of-f light constraint, the 

constraint function is found by adding the two powered-flight 

times, which are design variables, to the free-f light time as found 

from Equation A. 19, This value must then fall below the specified 

TOF (three years for most cases) for the constraint to be 
max 

satisfied. For the rendezvous requirement, recall from Chapter 2 
that the state vector at point i, x^, was defined such that 


j? 1 - if i 


(2.5) 



The rendezvous constraint requires that x^ = Xj TO , which 
represents the Emerald City ending the transfer with the desired 
velocity and radius vectors at the target orbit. The state vector, 
can be found from Equations A. 27 and A. 28, while the state 
vector, x^ can be determined from 


9 fro -e i + TOP (A. 29) 

N r P 



r x.iv - r p cos0 f _ ro 

^ y.TO “ TO 




(A. 30) 


v x,to ~ -v p sin 0 f<ro 
'V.ro " v- p cos0 f ro . 




(A. 33) 

With the 

above formulation, care 

must 

be taken 

in the 

programming of 

the above equations to 

assure 

that 

the 

proper 

quadrant is used when inverse trigonometric functions 

are 

used . 
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APPENDIX B 


ENGINE CHARACTERISTIC CALCULATIONS 
A computer program was written to calculate the important 
engine parameters for an open-cycle gas-core nuclear engine. This 
program is listed below. The eguations were either obtained from 
Ragsdale's paper or derived from basic considerations 


***************************************************************** 
program gascore 

***************************************************************** 

* this program finds the parameters of a gas core rocket for a 

* specific cavity diameter and Isp and prints the data to a file 

* called engine.dat. Enjoy. 

***************************************************************** 

dimension f(5) 

open ( uni t = 1 , f i 1 e= "ENGINE . DAT" ) 

GRAV=9. 80665 
PI=3. 14159265 
ATM=1 . 0132 5E5 
UGC=8314 
THRUST 3 4 . 4E+5 
VUVC= . 23 
write( 6,10) 

10 format ( '1' ,5X, ’ENTER SPECIFIC IMPULSE (Seconds)’) 
read* , SPIMP 
write( 6,40) 

40 format (5X, ' ENTER REACTOR CAVITY DIAMETER (Meters)') 
read* , DIAM 

** Calculate critical mass of Uranium 233 ** 


CD= 3 . 3*DI AM 

CRITM=201 . 836-67 . 0167*CD+7 . 70629*CD**2- . 158108*CD**3- 
& . 0196678*CD**4+ . 00088141*CD**5 

** Calculate mass flow rate ratio ** 

DOTM=5 . 041871E-5* ( VUVC**-10 . 799136) 

** Calculate hydrogen propellant temperature ** 

TH2=-294 . 461+4 . 36196*SPIMP+7 . 23731E-4*SPIMP**2-6 . 03953E- 
&7*SPIMP**3 

TH2=TH2+1 . 36261E-10*SPIMP**4-7 . 7 158 8E- 15* SPIMP** 5 
** Calculate reactor cavity pressure ** 
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PRES=( 1066. 57*CRITM*((THRUST*SPIMP)**. 277 )*(DOTM**. 1012)/ 
&(DIAM**3. 277) ) **1.383 


aa 


** 


aa 


aa 


AA 


AA 


AA 


AA 


Calculate mass flow rate from nozzle ** 

PDOT = THRU ST / ( SPIMP*GRAV ) 

Calculate mass flow rate of hydrogen and uranium ** 

UD0T= PDOT/ ( DOTM+ 1 ) 

H2D0T=PD0T-UD0T 

Calculate uranium temperature ** 

TTJ=140 . 06* ( ( PRES*THRUST*SPIMP/DIAM)** . 16 ) * (DOTM** . 017 ) 


Calculate pressure chamber volume and uranium volume ** 

VC=1.33333*PI*((DIAM/2)**3) 

VD=VUVC*VC 

Calculate hydrogen ratio of specific heats in cavity ** 

GAMMAC-3 . 84502- . 000336982*TH2+1 . 61689E-8*TH2**2-3 . 66763E- 
S13*TH2**3 

GAMMAC=GAMMAC+4 . 37181E-18*TH2**4-2 . 65532E-23*TH2**5+ 
&6.49491E-29*TH2**6 

Calculate hydrogen molecular mass in cavity ** 

H2MMC=2 . 32426- . 000111528*TH2+2 . 53066E-9*TH2**2-2 . 49052E- 
&14*TH2**3 

H 2MMC= H 2MMC + 8 . 94761E-20*TH2**4 
Calculate throat conditions ** 

TTH=TH2/ ( 1+ ( ( GAMMAC-1 ) /2 ) ) 

PTH=PRES/ ( ( 1+ ( GAMMAC-1 ) / 2 ) ** (GAMMAC/ (GAHMAC-1 ) ) ) 

PTH = PTH / ATM 

GAMMAT= 3 . 84502- . 000336982*TTH+1 . 61689E-8*TTH**2- 
&3.66763E-13*TTH**3 

GAMMAT=GAMMAT+4.37181E-18*TTH**4-2 . 65532E-23*TTH**5+ 
&6.49491E-29*TTH**6 

H2MMT=2 . 32426- . 000111528*TTH+2 . 53066E-9*TTH**2- 
&2 . 49052E-14*TTH**3 
H2MMT=H2MMT+8 . 947 61E-20*TTH**4 
ATH= PD0T*SQRT ( TH2 ) / PRES 

ATH=ATH/ SQRT (GAMMAT*H2MMT* ( 2/ ( GAMMAT+1 ) ) ** 

&( (GAMMAT+1)/ (GAMMAT-1) )/UGC) 

VTH= SQRT ( GAMMAT*UGC*TTH/H2MHT ) 

Calculate exit mach number ** 
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AMACH= 3 . 00 

50 ANSW=(2*(l+( (GAMMAC-1 ) / 2 ) *AMACH**2 )/ (GAMHAC+1 ) ) 

&** ( (GAMMAC+1 ) / (GAMMAC-1 ) ) 
ANSW=ANSW/AMACH**2-9E4 
if (ANSW.gt .0) goto 60 
AMACH=AMACH+ .05 
goto 50 

* Calculate exit conditions ** 


60 VE= SP IMP*GRAV 

TE= H 2MMC * ( VE / AHACH ) * * 2 / ( UGC*GAMMAC ) 

PE=PRES/ ( ( l+( ( GAMHAC-1) / 2 ) *AMACH**2 )**(GAMMAC/ 

& (GAMMAC-1 ) ) ) 

PE=PE/ATM 

AE=300*ATH 

GAMMAE=1 . 40018-3 . 60774E-6*TE-2 . 87154E-8*TE**2+5 . 90601E- 
&12*TE**3 

GAMMAE= GAMMAS- 2 . 94162E-16*TE**4 
H2MME= -2 . 8375E-6*TE+2 .01588 


** Calculate density of uranium ** 

RHOU=CRITM/VU 

** Calculate diameter of uranium ** 

DIAMU=2*(3*VU/4/PI)**( .3333333) 

** Calculate reactor power ** 

QR= ( ( TU/ 950 ) **6 . 25 ) *DI AMU/ PRES 
** Calculate thickness and mass of pressure shell ** 

THPS=PRES* ( DI AM+1 . 52 ) / ( 5 . 516E9 ) 

PSM=1333 . 33*PI* ( (DIAM+1 . 52+2*THPS)**3- (DIAM+1 . 52 ) **3 ) 

** Calculate mass of moderator, turbopump, nozzle, and 

** radiator ** 


EMM=191.67*PI*( (DIAM+1. 54) **3-DIAM**3) 

TPM= . 000307181*H2DOT*( ( 1 . 5*PRES)**( . 66667 ) ) 
ENM= 52861 . 26*THRUST/PRES 
PR=3E- 6*SPIMP+4 . 4E-2 
RM=145*PR*QR 


** 


Calculate total engine mass 


** 


ENGMASS=PSM+EMM+TPM+ENM+RM 


** 


Calculate engine specific mass and jet power ** 
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ALPHA 3 204 . 96*ENGMASS/ ( SPIMP*THRUST) 

PJET= . 5*THRUST*SPIMP*GRAV/1E6 

* Calculate exit velocity ** 

VEXIT=SPIMP*GRAV 

* Calculate thrust to weight ratio ** 

TWR= THRUST/ ( ENGMASS*GRAV ) 

PRES= PRES/ ATM 
write( 1,70) 

70 formatC SPECIFIED VALUES: ’ ) 

wri te( 1,80) SPIMP , THRUST , VUVC , DIAM 

80 format (IX, ' Isp= ' ,F8 . 2 , ' seconds ' ,3X, ’Thrust= ’ ,F11 . 2 , ’ 
SNewt ons ' , 3X , ' Vu/ Vc= ' , F3 . 2 , 3X , ' Dc= ' , F5 . 2 , ' Meters ’ ) 
write(l, 90) 

90 format ( // , IX , ’ CALCULATED ENGINE PARAMETERS : ' ) 
write(l,100) 

100 format(/, IX, 'ENGINE CAVITY CONDITIONS:') 
write( 1 , 110 )TH2 , TU 

110 format(lX, 'T(H2)=' ,F8.2, ’ Kelvin' ,19X, ’T(U233)= ' ,F9 . 2 , ' 
SKelvin' ) 

write(l , 120 ) PRES , CRITM 

120 format(lX, 'Cavity Pressure= ' ,F7 . 2 , ' Atm',13X, 

&' Critical Mass U233=’ ,F6.2, ' Kg') 
write( 1 , 130 )VC, VU 

130 format (IX, 'Cavity Volume= ' ,F5. 2 , ' Cubic 

&Meters' ,8X, 'Uranium 
&Volume= ’ , F5 . 2 , ' Cubic Meters') 
write( 1,140 ) DIAMU , RHOU 

140 format (IX, 'Uranium Diameter= ' ,F4 . 2 , ' Meters ' ,12X, 'Uranium 
&Density= ' ,F6 . 2 , ' Kg/Cubic Meter’) 
wr i t e ( 1 , 1 50 ) H2MMC , GAMMAC 

150 format (IX, 'Hydrogen Molecular 

&Mass( Cavity) 3 ' ,F5 . 3 , 2X, 'Gamma( Cavity )= ’ ,F5 .3) 
write( 1 , 160 ) 

160 format(/,lX, 'MASS FLOW RATES:’) 

write( 1 , 170 )DOTM, PDOT 

170 format (IX, 'H2-TO-U233 Mass Flow Ratio= ’ ,F7 . 2 , 6X, 

S'Propellant Mass Flow Rate= ' ,F6 . 2 , ' Kg/sec') 
write( 1 , 180 )H2DOT , UDOT 

180 format (IX , 'H2 Mass Flow Rate= ',F7.3,' Kg/sec', 6X, ’ U233 

&Mass Flow Rate= ' ,F6 . 3 , ' Kg/sec') 
write( 1 , 190 ) 

190 format (/ , IX, 'ENGINE DIMENSIONS:’) 

write(l , 200 )THPS , PSM 

200 format (IX ,' Pressure Shell Thickness= ' ,F5 . 3 , ' 

&Meters ', 3X ,' Pressure Shell Mass= ’ , f 8 . 2 , ’ Kg') 
write(l,210)EMM 

210 format(lX, 'Moderator Thickness 3 0.76 Meters' 

&, 8X , 'Moderator Mass 3 ' , F9 . 2 , ' Kg’) 
write( 1,220) TPM, ENM 
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220 format (IX, 'Turbo Pump Mass= ’ ,F7 . 2 , ' Kg' , 14X, 'Exhaust 

SNozzle Mass= ' , F7 . 2 , ' Kg') 
write(l,230) RM , ENGMASS 

230 format(lX, 'Radiator Mass= ' ,F9 . 2 , ' Kg' ,14X, 'Total Engine 
&Mass= ' ,F10 . 2 , ' Kg') 
write( 1 , 240 ) 

240 format ( / , IX , ' VARIOUS ENGINE PARAMETERS:') 

write(l , 250 )QR, .06*QR 

250 formatflX, 'Reactor Power= ' , F8 . 2 , ' MWatts ', 11X, 'Radiated 
&Power= ' ,F7 . 2 , ' MWatts') 
write(l ,260 ) ALPHA, P JET 

260 format (IX, ' Engine Specific Mass= ' ,F7 . 5 , ' Kg/KWatt', 

S3X , ' Jet Power= ' , F9 . 2 , ' MWatts') 
write(l,270 )TWR , P JET/ QR 

270 format (IX, 'Thrust to Weight Ratio= ' ,F7 . 5 , 10X, ' Engine 

&Ef f iciency= ' ,F5 . 4) 
write( 1 , 280 ) 

280 format (/ , IX, 'THROAT CONDITIONS:') 

write( 1,290 )TTH , PTH 

290 format (IX, 'Throat temperature= ' ,F8 . 2 , ' Kelvin' , 6X, 'Throat 

&pressure= ’ , F7 . 2 , ' Atm') 
write (1,2 95 )H2MMT,GAMMAT 

295 format (IX, 'Hydrogen Molecular Mass 

& ( Throat ) = ' , F5 . 3 , 2X , ' Gamma ( Throat ) = * , F5 . 3 ) 
write( 1 , 300 ) ATH , VTH 

300 format ( IX, ' Throat Area= ' ,F7 . 4, ' Sq Meters ', 11X, 'Throat 

&Velocity= ' , F9 . 2 , ' Meters/sec') 
write( 1 , 310 ) 

310 format (/ ,1X, 'EXIT CONDITIONS:') 

write( 1 , 320 )TE , PE 

320 format (IX, 'Exit temperature= ’ ,F8 . 2 , ' Kelvin' , 8X, 'Exit 

&pressure= ' , f 5 . 3 , ' Atm') 
write ( 1,325 )H2MME , GAMMAE 

325 format (IX, 'Hydrogen Molecular Mass(Exit )= ' , 

&F5 . 3 , 4X , ’Gamma(Exit)= ' ,F5.3) 
write(l,330 ) AE , VE 

330 format (IX, 'Exit Area= ' , F7 . 4 , ’ Sq Meters ', 13X ,' Exit 

&Velocity= ' , F9 . 2 , ' Meters/sec') 
write( 1,340 ) AMACH 

340 format (IX, ' Exit-To-Throat Area Ratio= 300 ', 10X ,' Exit Mach 

SNumber= ' , F5 . 2 ) 

* ************************* 

* This part of the program uses the method of Masser from Chapter 

* 5 to find the constant of integration for the calculation of 

* the dose obtained from the plume 

* 

f (1)= .07 
f(2)= .00002 
do i=l , 2 

BIG=GAMMAE*SQRT( 2/ ( GAMMAE- 1 ) * ( 2/ ( GAMMAE+1 ) ) ** 

& ( (GAMMAE+1 ) / ( GAMMAE- 1 ) ) ) 

CFCFMAX= ( BIG*SQRT( 1- ( PE/ PRES ) ** ( ( GAMMAE- 1 ) /GAMMAE )) + 
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350 

360 


&300*PE/PRES)/BIG 
X=l/ ( SQRT( PI ) * ( 1-CFCFMAX) ) 

B=X/4/SQRT(PI)*( (GAMMAE-1)/ (GAMHAE+1 ) )** . 5*(2/ (GAMMAE+ 
Sl))**(l/ (GAMMAE-1) ) 

TOP=f (i)*QR*6. 448E+7*(H2MME*PE/TE/82.06)*AMACH*B* 
&AE/PI*10000 

BOTTOM 3 H 2 DOT *3 . 7* ( 1+ ( GAMMAE-1 ) /2*AMACH**2 ) ** . 5*( 2/ 
&(GAMMAE+1) )**( (GAMMAE+1)/ 2/ (GAMMAE-1)) 

CON ST= TOP / BOTTOM 

format (//IX, ' lambda 3 ' ,F15 . 9 , 5x, ' constant 3 ' , f 20 . 4) 

wri te(l, 350 )X, CONST 

format (IX, ' f ( t ) 3 ' , f 10 . 7) 

write(l,360)f (i) 

enddo 

close(unit=l) 

stop 

end 
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30 ssfiBqusoJsd iBnqoB sqq ' SAoqB psqxo suosbsj owq sqq 203 sqsuBjd 
2suux sqq oq sasqsusaq sqq 30 sanqBU sqq 30 ssnEOsq ' snqj, 

•qi^ssa 2sqqsq qonui 

b 's/uo[ 8 ‘ 81 sq A juo pxno« AV ©qq 'S2BS R £fz ^Tuo ux JisqsuBaq sqq 
buxqBUi Aq 'asASttOH *s/un( sq pinow psaxnbsa AV sqq ' S2BS/ 5*3 

30 suixq qqbxx3 b qqx« 2S3SUB2q b S2{bui oq /2q oq sjsm Aqxo px*.3©ui3 
sqq 33 'quxod sxqq qv 'S2BSjt 5*3 uxqqxw S2sqq sq oq ps2sp2o 
sx puB 'OSOZ ‘ LZ itJBnaqsj uo snusA oq pssoojtd oq saspjo ssaxsosj 
Zqxo px*-*sui3 sqj, 'dWI'inW ^q p©qoxps2d ' sxduiEXS fiuxMonoq sqq ux 
usss sq ubo uousuiousqd sxqj, • sn x ba psxqxosds sqq sq oq suixq qq£fX3 
sqq 203 2sp20 ux sqnoa qosaxpux ue buxipeq £q ix^s IIT M 3 n q 'qsuBXd 
qsbJBq sqq oq Ax 3 D9:i TP 0)5 3 ou IIT« AJoqosCBaq sqq ssriBDsq sx sxqj, 
* sssbsjdux Axx^qoB psjxnbsj AV ©qq ' S2BSA ss2qq ' Abs 'ux sqsusxd 
sssqq 30 suo oq dX2q b fiuxojoq £q 'snqj, -sjbsA o«q oq suo 30 
jspjo sqq uo sx sqsusxd 2suux sqq 203 (AV qsswox ©qq UT buxqxnssj 
suixq qqfix x 3 ©qq) ©uixq qqbxx3 x BUJt T 3 <*° sqj, -sqsusxd 2suux sqq 
203 ssbo sqq qou sx sxqq 'asASMoq ' x^asusb ui -uoxqBjnp fiuxpuBUisp 
sssx ©qq 30 ssnEOsq ' 2 sbuox spsui sx suixq qqbxx 3 ©qq sb u«op 
06 pxnort AV sxqq qBqq qosdxs qqfixui suo ‘punoq sx JsqsuBjq qoBxs 
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•Aq pejmstsui ©q ubo sqxqjo jbutuiou quejejjxp 
qcqq spxqspjtA e ' qutqjioduix seuioosq AV P sx T3 ® Jog ©xqxssod si 
jegsutxq t qtqq auixq 50 aPtqueojad ©qq 'sqautxd aeqno ©qq Jog ' snqj, 
•suoxqxpuoo umuixqdo ueqq ss©t qt saegsutJiq ©3 [bui oq pexedead ©q qsnui 
put 's©t6ut ©stqd © xqxssod ©qq 50 uoxqotag t Ajuo j©quno©u© XIT M 


U'TT 

00’ 3T 

16 ’31 

30 ’ 81 

99 ’961 

nv s = « 

LI' 8 

Tfr ’8 

£8 ’ 8 

86 ’0T 

39 '\Z 

nv > = « 

ic’s 

L£ m S 

frS’S 

T£ ’ 9 

93*6 

nv z = u 

(sjb©A) uis i 

oqnxd 

aunqdeu 

smiBjfl 

xun^BS 

jaqxdnf 

qautXd 


sqeut xd Jeqno eqq 5 ° spoxjad oxpouAs - ZT'Z 
Aqxo pxtJsuia aqq 'sxqq 50 asntoag "ZT’Z eiq^J. ux uaas st 'sqsutxd 
jsqno ©qq jog a© 6 uox qonui ©xt spoxjad oxpouAs ©qq asntoaq sqautxd 
jsuux ©qq oq utqq Axquajaggxp peq«sJ3 ®q qsnui sq©u©x d 3 ©qno ©qq oq 
saegsutaq ©qq go uoxqt 6 xqs©AUX ©qj, *s,n¥ 01 P UB 'S ‘f '£ 3° sqxqao 
X tuxuiou jog OT’Z-L'Z ssxqTBj, ui punog ©q ubd (oqnxd oq jaqxdnf) 
sqsutxd jaqno ©qq jog sxsAxeut uoxssxui ©qq uiojg sqxnsaj ©qj, 

sqautxd Jsqno 

♦ 

•sjtaddt uoxqtquexjo ©xqsjOAtg ©qq XT3^n qxtM uto qx 
asntoaq 'sjapjo sqx 6uxax©o©j uiojg sjtaA ©ejqq japun ux auixq ©qq 50 
% 00 T spsuJ aq uto jagsutjq aqq qtqq pxts aq uto qx uaqq 'Ay u©ax6 b 
qqxM ss©x Jo sjb©A 98 ‘T U T ©uixq aqq 50 quaojad Aut paqotaj aq uto 
qqjta gx 'snqi -poxjad oxpouAs aqq 50 qqfiuax aqq oq ©np psjjnooo 
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paqoaax aq ubo ujnqas •AxsnoxAa.xd sqauBXd aauux aqq 30 uoxssnosip 
aqq ut dn qqfinojq ajaw qaqq squxod aqq Puxxapxsuoo 'ssax xo sJBaA 
aaaqq ux s/unj OV 3 ° S ,AV paqoBa.3 aq ubo sjbw pub qq.JE3 ' snuaA 

'osxv -AqxxxqBdBO AV tpnui sxq^ qssax qa apxAoxd oq aAaq pjno« 
uiaqsAs uoxsxndojd aqq ' qaui aq oq sjb jaqdBqo sxqq 30 fiuxuuxfiaq aqq 
qa u«op pxB x sxsn qaqq sxeofc aqq jx ' snqj, ’S/uni ofr jo AV ^ B «- auo 
a qqxw auixq aqq 50 % 00 T paqsx xduioooa aq ubo aaqxdnf oq aajsuajtq 
jaaA-aaaqq b ' qxqjo sxqq uioxa -jaqxdnr pua sxbw 'Axnquao ^XZ 
-pxui aqq ux AqxAxqoB aoads AxaqauB xdxaqux 30 squxod x boo 3 paumsaad 
aqq oq asoxo pua 'punoq aq uao sx'ex.xaqeui axqnnxaA 30 sqxsodap 
qsBA ajtaqw 'qxaq pxoxaqsa aqq oq jaau Ax^AxqaxaJ Aqx:) px^-xauia aqq 
sqnd qxqao sxqj, -nv 2 ’S 3 * 'aaqxdnf puB flV S'T 3 ® 'sjbw uaawqaq 
qxqjo sxqq 30 uoxqeoox x 12 - 3 ^ 11 ® 3 sx ® uo 3 s:i T 3 sqj. 'saoqoB} 

XBjaAas uo pasaq sen uoxqoaxas sxqj, -qxqao x BU T mou Aaauxuix xaJd 
aqq sb paqoaxas sew s,nv 3 ® qxqao UE 'uoxssnosxp auios puB dWHQH 
Aq paqaaauafi saxq^q afiaquaoaad aqq J 5 uxuxuiaxa aaqqv 'S.nv OT pu« 
' g 'fr '£ a jaw pauxuiaxa ajiaw qaqq sqxqjo x^T 1110 * 1 ©X dulES 

uoxqoaxas qT^O X BU T UIO N fr' 1'3 

• qxqxo 

Xbuxvuou aqq 50 uoxqoaxas aqq X03 uoxqajapxsuoo quaqaoduix ub sew 
sxqj, "uwop quaH Axx Bn 3 DB saSaquaoxad aqq ' saseo auios ux puB 'uns 
aqq uiojj Aewa qxqjo x BU T m0U ®*n PuT AOU1 pauxafi qxjauaq axqqTI 
sbm axaqq 'OT 'Z-L'Z s®I<T B l ui uaas sv "sqauaxd xaqno aqq qoaaj 
oq Aqxo px^Jauia aqq 30 Aqxxxqa aqq uo paq uns aqq uio.33 aaqqanq 
qxqjo x BU T UIOU a*n Puxaoui qoaqqa qaqw paqou aq pxnoqs qi 
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suit:* sqg pus 'sxsdBttsd sqg 10 gusumbtB sqg 'spou buxpusosB sqg go 
spngxbuox sqg 'uoxgBuxxoux sqg 'sxxe toCBUi-xuiss sqg ' A*iou;uaooe 
aqg .'sgusuisx© x^^T^- 30 Ihdtssbxo xxs sqg buxsn Aq pstsgus sx Agxo 
pXEtsuia sqg go gxqto sqg 'spoo dWimw sqg buxsn usqM -uoxgdumssB 
sxqg go AgxpxxEA sqg og sb psstBt sas« sgqnop suios ' uoxgBbxgssAUt 
tsqgtng uodn 'xsasmoh -AtBSSsosu sq pxnoM gBqg AV umuixxBui sqg 
©Xqnop ^x 3 ! 6310 - 3 PI^om dxtg utngst sqg tog buxgunoooB gBqg psumssB 
AXIEtousb sbm gx 'Apngs sxqg ux Ax^^a -gxqto X® 33 ™ 1031 s3 ig °3 utngst 
sqg sx sptAOtd gsnui uisgs^s uoxsxndotd sqg gBqg Ay AtBsssosu sqg 
buxssnosxp usqM dn gqbnotq sq pxnoqs gBqg guxod buxgsstsgux uv 

•(sxqxbgxbsu Zxx® 3 ® 31 ® 6 
stsM gno gxqto sqg Buxaoui £q psutBb sgxgsusq sqg) psgosxss gou sb m 
gxqto snxpBt tsbtBX e Aq« sx qoxqM ' ssxousbtsuis puB suoxgBoxununuoo 
tog qgtsa go sousgsxp sxqBUOSBSt uxqgtM Agxo pxEtsuia sqg sgnd osxb 
s.jqv ?o gxqto xeuxuiou V * ssxgxxTCfEdBD uoxsxndotd sqg buxssstoux 
to stsqg Abm sqg go gtBd buxxsngst Moqsuios sb qons ' suoxgBtspxsuoo 
XBtosds stxnbst pxnoM sgsusxd tsgno sssqg og sdxtg ' snqj, -s/uiq ofr 
go gxuixx AV e pub stssA sAtg og dn ssuixg gqbxxg qgxM stsgsuBtg Zq 
X X e gB psqoBst sq gouuBD ognxd puE sungdsw gnq 'pstxnbst sx gsuBXd 
gBqg og dxtg b gx suixg sqg go %00T psqosst sq ubd Atnotsw -Atnguso 
qj-pxui sqg ux uoxgBtoxdxs tog guBgtoduix AxXEUtbtBUi psuissp sgsuBxd 
go ZtobsgBO sqg ogux XT®3 snuBtn ©snBOsq 'sxqEgdsooB sq og gx©3 SBM 
sxqj, -suixg sqg go %££ Ax uo ©X^T 5503 S T -isgsuBtg sqg ' stssA ©Axg 
tog 'puB 'pssn sq gsnui stBsZ SAtg puB tnog go ssuixg gqbxxg 's/uiq 

Ofr go AV e qgxM snuBtn qoESt oj, • snuoq psppB ub sx qoxqM ' stBsA 

sstqg go ssuixg gqbxxg buxsn s/uiq ofr 3 ® suixg sqg go %00T Ax-iesu 
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ux snuEan qoBaa oq %os aAoqB pus sjbsA aaxqq ux uxnqBS puB jaqxdnf 
qoBaj oq AqxxxqB %00T ux qxnsaj Pl^oqs aanbxg sxqq qBqq qx©3 
sen qx 'sqauBX<i Jaqno aqq joa * aAoqB passnosxp sbm sb saagsuBjq 
XBUixqdo sjfBui oq AqxxxqB s ( Aqxo piBjauia aqq oq anp eqsnbapB aqxnb 
aq pjnoqs axnbxg sxqq 'sqauBxd xsuui aqq aoa -uiaqsAs uoxsxndojd 
aqq oq aAXb oq qxuixx jaddn AJBSsaoau aqq sb pauinssB sbm s /uiij 
OS 3° AV I E 3°3 u1s os 'ST^q S3 (TI uoxqBnqxs b jog AV ©qq aqs^noxBo 
oq pasn aq pxnoo dHIinw woq jo gx u«oui[ qou sbm qi -paqoadxa 
A X TBUXbxjo sb qbxq sb aq qou Axqsqojd pxnow s,AV ©qq puB 'jagsuEjq 
aqq jog AV qsBax aqq ajxnbaj pxnon qsqq quxod aqq qB qxqjo x^uxuiou 
aqq oq uanqaj pxnoo Aqxo px^Jauia aqq 'uiopaajg sxqq go asnsoag 
•qxqao qsqq uo axaqttAuB go psaqsux uoxqBaox M s t quxod qsoqb u 
aqq oq ujnqaj oq psq qx gx sb qx aqBxnox E£) Pinow qx 'dxjq ujnqaj 
aqq aog pajxnbaj Ay E qno ajnbxg oq pasn aaaM awnnw 31 - qxq;io 
qBqq ux qauexd b aaaw qx gx sb 'qxqjo x BU T ulou peuxgap ©M3 punojB 
saAoui qBqq paqaejo sx 4l quxod qsoqb,, b 'qxqjo x^uxuiou ©*13 ©uxgap oq 
pasn ajB squauiaxa x^qT^- 30 ©^3 uaqM ' Ax©qBunqaogun "qxqjo x eu T ulou 
aqq 6 uoxb quxod Aub oq ujnqaj oq aaag sx qx 'psaqsui *qxqjo 
XEUxuiou aqq uo quxod JBxnoxqjBd b oq uanqaa oq SBq Aqxo px®- 3 ©^ ©qq 
sAbs qsqq buxqqou sx ajaqq 'dxjq uanqaj aqq jog 'X9A9MOH 'dWIinW 
Aq paqBxnoxBD aq ubd AV ©qq puB 'qauBX<3 qabJBq aqq qB baxjjb 
pus qxqxo x® 313 * 11011 ©qq uo quxod JBxnoxqjBd b uiojg aABax IIT W Aqxo 
pXBjauia aqq 'qauBxd qabJBq b oq qxqjo x EU T UIOU ©M3 uiojg buxob uaqM 
*auixq uaAXb b qB qxqjo aqq ux Aqxo pxBjauia ©qq saqBoox qsqq enxBA 
auios Aq paoBXdaj sx quauiax© qsBX aqq puB 'sxsdBXjad ou sx aaaqq 
'pauxuiBxa aaa« qsqq sqxqjo JBxnojxo aqq joa *abBSSBd sxsdBXjad go 
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• *oqa 'spoqqaui aouajagggp-aqgugg Aq ao Ax X BO T}^I BUB paqanxcAa aq oq 
aaE suogqoung aqq go squagpajb aqq jaqqaqw 'jaaugx-uou jo jaaugx sg 
uiaxqojd aqq aaqqaqw ' squgajqsuoo go adAq aqq iuiaxqojd eqq go ainqeu 
oqq uodn spuadap paqoax©s poqqaui jaxnogqjad aqj, -uuog 9 AoqE aqq ug 
suiaxqojd uogqazguigqdo 9 ajos oq sanbguqoaq uo padoxaAap Aaoaqq qonui 
xiaaq saq ajaqj, ‘©^sq ued saxqtxjEA asaqq qaqq sanx©A go abuaa aqq 
sagggoads puB saAxasuiaqq saxqsgJBA ubgsap aqq uo sg qugajqsuoo go 
adAq pjgqq ©qj, * squgajqsuoo Aqgx^nbaug paxx® 3 S T ©<*^3 puooas aqq 
axgq« squgaaqsuoo Aqgx©nba paxx BD feuxaq adAq qsagg aqq ' squgajqsuoo 
ui ajE aaaqj, • sa xq^gJUA ubgsap u *a'X 'u qqbuax go joqoaA axqsgjBA 
ubgsap aqq sx x puB uogqoung aAgqoaCqo aqq sx (x)g uogqoung aqi 

• n x > X > lx 

w' • • • • "C+ 3 U i = C jog o < ( x )^ 6 
9 ui' ,, - ,/ x = C jog o = ( x )^ 6 

squgajqsuoo aqq oq qoaCqns 

(3-3) (*)J '* 9X uxui 

Aq paquasajdaj sx uiaxqojd sxqj, •saxqug.ieA sqg 
uo squgajqsuoo ugaqjao oq qoaCqns 'ugaqqo ued 'uogqoung aAgqoaCqo 
aqq paxx B3 'uogqoung a umuigxaui ao umuigugui aqq buguguuaqap go 
auo sg uiaxqojd uogqazguigqdo ogsaq aqj, *ssaoojd ubgsap bugjaaugbua 
aqq ug x°°l x n 3- 3 ® wod B S T uogqazguigqdo go PX®T5 ©q.1 

uogqonpojqui XTZ 
NOIiLVZIWINIW AV Z'Z 

• suogqajapgsuoo x® 30 ®* 35 ©uios qnoqqgw axqsqoaajun aja XIT3 S 
oqnxd pua aunqdaN 'asjnoo go ’(sasao auios ug 'ssax jo) sjaaA aAgg 
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buxuuogaad *p©sn ©q pjnoD qBqq s©x6ub qsnaqq ©qq ao 'pejxnbsa 

qsnaqq 'sauixq qqbx jg-peaertod ©qq aog pjBbea ou qqx« 'ZqxooxsA 
ui ©bUBqo snoauB^UE^sux ue pauinssB qBqq uiBJboad asqnduioo e pasn 
aaeqq psquasead sxs^xbub ©qq 'jsasmoh 'aeqdBqo sxqq ux aaxx-JB© 
pequesead sbm q©qq qjow ©qq qqxw quEpunpea ©q pxnow sxqq qBqq uiaas 
i^Bui qi ’©xqxssod sb ©xqBuosBea sb squauieaxnbsa uieqsits uoxsxndoad 
da©3j oq aapao ut pexgxqsnL sb« AV ?o uoxqBzxuixqdo ub qBqq qx©3 sem 
qx ' uo sb© j sxqq aoa ’©xqxssod sb wox sb s,AV ®s©qq d©aj[ oq Aqx q 
piBJsuia ©qq go qssasqux qseq ©qq ut sx qx ' snqj, ’'-oq© 'abEaoqs 
quBXIBdoJd ' so xuiBuAp x 12 - 311 '} 013 - 3 } 5 'buxpxsxqs ut sui©xqoad qxnoxggxp 
©aoui oq pB© x ( 1 * ‘sauixq qqbxxg-psaaMod 'sqsnaqq) squsuisaxnbaa uuaqs^s 
uoxsxndoad aeqBaab Aqxo px^J^uia ©qq go suisqsAs ©qq go qsea ©qq uo 
qoegg© apEosso b SBq sxqq 'uieqsAs uoxsxndoad ©qq buxuxnaqs sapxsag 
• ©XqBjapxsuoo ©uiooaq uiaqsAs uoxsxndoad ©qq uo spuBtusp ©qq '©SBsaoux 
AV aog senxBA ©qq sv * squeuiaaxnbsa uiaqsAs uoxsxndoad ©qq 

uo qoEduix qB©ab b © ABq saegsuBaq ©saqq aog p©p©au AV e q3 q^qq psqou 
ua©q ssq qi *ui©qsAs aBXos ©qq ut sqsuBXd ©qq oq qxqao x^™ 011 B 
uioag saegsusaq ©qq qqxw sxB©p HSIM qostoaa go qusuidoxsAsp go saBsA 
ortq qsaxg ©qq ut peoBgans SBq qBqq sui©xqoad aoCBui ©qq go ©uo 

punoa^3{OBa ZTITZ 

• uoxqBZXuixuxui 

A V 5° uiexqoad ©qq ^x 9UI1BU 'pessnosxp ©q XIT W JB©A sxqq pexpnqs 
ui©xqoad uoxqBzxuixqdo ub 'puxui ut uoxssnosxp © AoqB ©qq qqxM 

•uiexqoad uoxqBzxuixqdo ©qq ©axos oq pasn ©q oq ©poo asqnduioo 
©qq 6uxqo©XBS uaqw paaepxsuoo ©q oq psau XIT M saoqoEg ©saqq go x IV 
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'©sjnoo esfr ^-©-[T^O 0 Coad b «on°? og psumssB sx Agxo piBasuig ©qg 
'ssuixg gqbgx g-ps-iewod ©qg PutJna ‘{ Z'Z ©Jnfexg ©es) AgxAB.i£ s t uns 
©qg go lind ©qg Aq pegoeggsun st seuixg gqbx xg-p©.J©Mod ©qg 6ux;xnp 
AaogosCBjg ©qg gBqg sx uoxgdumssB ' x©goAXd gsoui pus 'pxxqg ©qj, 

•©SBqd gqfex xg-ps-JSMOd puooss ©qg 
go 6uxuux6©q ©qg gB uoxgBJ© x©oo© ©qg og x Bn ^® g©s sbm ©SBqd gqbxxj 
-p©j©«od gsjxg ©qg go pu© ©qg gB uoxgB.x©x©oo© ©qg 'uoxgdumssB sxqg 
6uxx[bui Ag -Apngs AxBUXUixxsJd sxqg aog ©xqeoxxddB ©q XIT^ S PinoM 
sgxnsej ©qg gBqg peuosBSJ sbw gx gnq ' © xq©uoxgs©nI) sx uoxgdumssB 
sxqg go Agxpxx©A ©qj, -stsAxbub sxqg xog pexouSx sbw gx ' Agxoxxduixs 
jog gnq ' squBg ©bBJogs ©qg qbnojqg ©6 bj[b©x gB©q og ©np usbojpAq 
pxnbxx ©qg go ggo-xxoq ©uios ©q XIT M ©-asqg '©sanoo go *a©gsuBag 
©qg go ©ssqd gqbxxj-ssaj ©qg buxjnp guBXX©doad go ssox ou sx ©asqg 
gBqg sb« ©psui uoxgdumssB puooes ©qj, 'HSIM goeCoag ux gssasgux 
gsB©x ©qg jo ©jb qoxqw 'ognxg puB Aanojsw gdeox© sgsusxd ©qg go 
XXB aog uoxgdumssB ©xqeuosBea b sx sxqg ‘ T'Z ©iq©i ux Urtoqs sbm sv 
•qgjBg ©qg sb ©UB]d ©uibs ©qg *©"x '©UBXd oxgdxxo© ©qg ux xt B ® J1B 
puB sgxqxo jbxhojxo ux xi® ©2© sgeuBXd gBqg sx uoxgdumssB gsjxg ©qj, 
* ©psui ©q og psq s^egsuBJig puB sgxqjo ©qg go ©jngBU ©qg fiuxpjBbsj 
suoxgdumssB uxBgxeo 'uoxgBzxuixgdo ©qg gno Axjbo og Jtepao ui 

suoxgdumssv £’3*2 

*AV gs©Mox ©qg px©T^ IIT W 3©qg 
'gsuBxd uoxgBUxgsep ©qg jb©u Axx© 0- ?* 1 *} 'JTq- 10 g©6J©g ©qg pus Agxo 
pX©J©uia ©qg u©©«g©q uoxgBgu©xjo ©qg buxuxuusgsp sb xi® w S1B 'gunoooB 
ogux uaqBg ©q ubo sJogoBg ©seqg gBqg pedoq sx gx 'uoxgBzxuixgdo sxqg 
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•UOXJB6XJS©AUX AJBUXUIX X 9 J d SXqj JO} PTT EA SBM UOXjdUmSSB STqj JBqj 

1 XSJ sb« j x snqj pus 'jequmu ©6 :ibx Ajsa b ©q xiT tt h '^ 3 TO piBjauia 
©qi joj ssnxEA 8AT}E}uessJdej JOd * (3 )ui ' ssbui dxqs puB '(Jp 
'snxpBJ ©qj ' 1 \g 'jsnjiqj ©qj uodn juapuadap sx h 3®q3 usss ©q ubo 31 


(C’ 2 ) 


— XT 0I X LZf I- Tl 

(3 )arri f g 
t (3 ) " 173 P 


21P-°UI m 

■ ( 3 ) 2 T 


( 3 )W 


Aq u©ax 6 'sessqd jqfixx J-peJSHOd ©qq fiuxanp 
AqxABJfe puB jsnaqj 03 ©np suox^bj© 1 ©oob jo oxjbji ©qq jo ©r^ba 
© qq uo pasBq sx uoxjdumssB sxqj, • sjiojo©a (qsnaqj) uoxjbji© xsoo® 
puB Ajxdox©a x«T3Tux ®q3 3° sjuauoduioo-A puB -x ®qq peuxuuejsp 
sasBqd jqfexxd-peJ9Moa -xoj goxqdumssy - £’g ajnfexj 
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qa qTqao qq6tm-esJ5 30 
qa qTqjo qqbTxq-aa^q 50 


3 quTod 
£X bu, oub anaq 
X quxod 
£X bui oue anjq 
(sapxsdB jo 
auxx ) aoqoaA AqxoTjquaooa puB stxb-x 
auieaj sousisjsj X B T^ JSU T uaewqaq axbuB 
qounsx qa qauaxd qabJeq puB stxb-x 
auiExq 93U9J9)9J x B Tq j ® u T uaaMqaq ©X^ue 
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